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ABSTRACT 


This is the second volume in a two-volume final report series for 
Project NAS 9-14467 sponsored by the Earth Observations Division, NASA/JSC. 

This report series summarizes the work covered between the period November 
15, 1974, and November 14, 1975. The objectives of the project were to 
evaluate the LAC IE II table look-up approach to sun-angls correction. Canopy 
reflectance modeling was employed as a technique for evaluating sun-angle 
signature extension. 

Volume I presents the multiplicative and additive coefficient matrices 
for a linear sun-angle correction approach. These coefficient tables are 
calculated using either measured empirical canopy reflectance functions or 
model derived data. These values are then incorporated into an atmospheric 
radiation transfer model. The dependence of the coefficient matrices on 
crop stage, crop type, and canopy directional reflectance variations is 

reviewed. Finally, a method for inferring leaf area index, an intrinsic 

■» ► 

scene characteristic, from canopy reflectance is discussed. 

Volume II presents the basic data and computer programs used in the 
study. A brief review of the radiometric and geometric data collection 
procedures is also given. In particular, two recent methods developed by 
the investigators for determining plant geometry are discussed. These include 
the Fourier diffraction and multiple view angle approach. The data compila- 
tion consists of canopy reflectance, constituent reflectance. Leaf -Area- 
Indices, and leaf slope distributions for four wheat crop development 
stages at Garden City, Kansas. 
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I. INTRODUCTION 

This is the second volume in a two-volume report series for Project 
NAS9-14467 which represents Task 4.1.1.2f(s) in the LACIE 00200, Volume VIII. 
Specific objectives of this task in order of priority include: 

A. To evaluate the current LACIE II table look-up sun angle correction 
algorithms relative to; 

1. The effect of canopy reflectance variations with sun angle; 

2. The effect of canopy sun angle reflectance variations with crop 
development stages; 

3. The effect of applying a uniform sun angle correction developed 
specifically for wheat to all crop types. 

B. To recommenc modifications to the current LACIE II sun angle correction 
algorithm. 

C. To investigate the alternative sun angle correction procedures for present 
and future satellite systems. In particular, to investigate the 
possibility of extracting intrinsic scene characteristics from wheat 
canopy modeling. 

The first volume in this report series summarizes the alpha, beta coef- 
ficient matrices required for a linear sun angle correction algorithm. The 
effects of crop stage, type, and canopy directional reflectance properties 
on the correction approach are reviewed. In addition, the relationsh-ip of 
Leaf-Area Index to canopy reflectance is discussed. 

This volume provides a brief description of the field data collection 
techniques. Section II, and a complete listing of the radiometric and geometric 
data measured at the Garden City, Kansas, test sites for each of the four crop 
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development stages* Section III. Program listings and control cards for the 
computer programs used in the study are given in Section IV. The recently 
developed Fourier diffraction and multiple-view angle techniques for assessing 
plant geometry discussed in Section II should be of particular interest. 
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. FIELD DATA COLLECTION TECHNIQUES 
1.0 Radiometric Measurements 

All of the reflectance data given in Volume II and analyzed in this report 
were collected with an Exotech ERTS radiometer which was modified for digital 
read out. The field procedure involved taking two readings centered on a row, 
then two more centered between rows for a total of four tripod movements. 

A wh-i .e panel reference reading was made at each set-up. A final measurement 
of diffuse/direct irradiation was made by reading a sun-lit and shaded 
reference panel, respectively. These measurements were repeated for each of 
three intensive plots four times during the day. In addition, periodic bare 
soil reflectance measurements were made. Figures 1, 2, and 3 show the 
radiometer in its various field configurations. 

Leaf transmission measurements were made using a special attachment for 
the ERTS radiometer that was developed by the Laboratory for Applications 
of Remote Sensing at Purdue University. The attachment consists of a 
cylindrical barrel, a sphere, and two flat discs which have small slots cut 
in them between which the leaf or other material is placed. The apparatus 
fits over the individual sensor ports on the ERTS radiometer and must be 
shifted to measure transmission for each of the four bands. Once the radiometer 
and attachment are aligned with the sun, it is a simple matter to shift from 
port to port for measurements in each band. The unit is aligned fay the use 
of a pin-hole type site on the side of the barrel assembly. The barrel is 
directed towards the sun and collects the direct rays which pass into the 
sphere. The sphere is coated on the inside with BaSO^ and contains a 
blocking baffle which prevents the passage of direct solar radiation through 
the sphere, allowing only diffuse radiation to pass. The diffuse light then 
passes through the material contained between the disc holders and is recorded 
by the sensor. The attachment is illustrated in Figures 4 and 5. 


original page lb 
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Fiaure 1. ERTS radiometer 
in' its field configuration. 



Figure 2. Four measurements were made 
(2 over a row and 2 between rows) for 
each "set-up" at an intensive site. 

Four such set-ups were completed during 
the day, covering a wide range of sun 
angles. 
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Figure 3. Percent reflectance for 
each band was determined by dividing 
the sample reading by a reading made 
from a white reflectance panel, which 
approximates total scene irradiance. 




Figure 5 • A baffle in the center of the sphere prevents direct sunlight from 
hitting the radiometer's detector and saturating the signal. 

PAGE IS 
OF POOR QUALnyi 


Figure 4. The leaf transmission attachment, in turn, is placed over each ap»''^ture 
Two readings are made; the first is an unobstructed measure to establish a 
reference signal with the second having the leaf placed over the apeture of the 
attachment. The percent transmission of the leaf is calculated by dividing the 
reference signal by the sample reading. 
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All radiometric data were recorded manually in the field as illustrated 


on the enclosed data sheet. 
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2.0 Techniques for Assessing Leaf Angle Distributions 

A. Introduction 

The Second quarterly Progress Report referenced our preliminary work in 
developing rapid in situ techniques for measuring leaf slope distributions 
in wheat canopies. In particular, we discussed the Fouriev’ diffraction and 
Fredholm approach. The field procedure for these two techniques has been 
implemented as part of the LACIE Field Measurements Program at Garden City, 
Kansas. Two alternative methods for measuring leaf slope distributions include 
the standard point quadrat approach and the orthogonal photographic method. 

The procedures for all four methods are described below. 

B. Point Quadrat Technique 

It has been shown that the mean foliage angle can be calculated from the 
number of contacts made by point quadrats passed vertically and horizontally 
through a plant canopy (Wilson, 1959). In practice the error associated with 
this method rarely exceeds 10?5 (Wilson, 1962). The technique is in situ , 
however appreciable localized trampling is induced around the field plot. The 
time required for a single angle determination is about 18 man hours (Knight, 
1973). This method is most commonly used to characterize foliage geometry, 
however, it only estimates the mean inclination angle rather than a 
distribution. 

The field procedure involves the calculation of the average number of 
contacts a long slender pin makes with the vegetation during a pass through 
the plant canopy. The length of the pass, for both the horizontal and vertical 
transects, is dictated by the height of the canopy. Several hundred passes 
are made from both directions, with the averages for each being multiplied by 
theoretically obtained coefficients to determine the mean foliage angle. 
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C. Orthogonal Tracing Technique 

The distribution of foliage angles for an individual plant can be accurately 
determined by analyzing two orthogonal photos of the plant (Oliver and Smith, 
1974). The distribution of angles for an entire plot is statistically 
determined by averaging the distributions of several representative plants. 

This technique has many of the limiting features associated with the point 
quadrat method. It is slow, tedious and destructive. However, its accuracy 
makes it a prime technique for evaluating the results of the other methods. 

With this procedure, individual plants are clipped from a field plot 
and the silhouetted profiles are photographed from two orthogonal directions 
(Figure 1). The photographs (Figure 2) are then digitized by placing a 
transparent grid over the photographs and recording the two-dimensional 
coordinates of straight line segments along the profiles. The profiles 
are then plotted on microfilm (Figure 3) using the digitized data in 
order to verify the hand digitization. A computer program determines the 
three-dimensional coordinates of the foliage elements from the two sets of 
orthogonal data, and calculates the average foliage inclination angle by 
direct computation (Figure 4). The distribution for the entire plot is 
calculated by weighted averaging of the individual plant distributions based 
on the size of the plant. 

D. Fredholm Integral Technique 

This technique has the easiest data collection and reduction procedures 
of all the proposed methods. However, the technique cannot be applied to 
dense plane canopies (Oliver and Smith, 1974). 

The proportion of gap or "probability of hit," as a function of view 
angle is a function of the mean canopy projection in the direction of view 
averaged over all foliage elements. Several explicit expressions of this 
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Figure 1. The field procedure for the 
orthogonal method involves making 
orthogonal photos of a silhouetted plant 
Markers for several branches are used to 
avoid confusion when a pair of photos is 
digitized. 




Figure 2. Orthogonal 
views of a silhouetted 
plant. 




Figure 3 . Computer 
plot of the digi- 
tized photos. 
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Figure 4. Distribution 
of angles of the plant in 
three-space. An estimate 
of angular bias within the 
canopy is made by averaging 
the distribution of several 
plants . 
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functional dependence are given in the literature (Nilson, 1971). For 
example; ®r 

where' g(9r) is the mean canopy projection in direction 9^.. LAI is leaf area 
index. 

Given a measured P^tOr), we can then invert the expression to derive 
g(9r). 

This mean canopy projection in direction 9r, g(9r) can then be related 

to the leaf slope distribution f(0a) via a Fredholm integral equation of the 

ir/2 


first kind (Oliver and Smith, 1974) g(9r) 


' I 

J 0 


K(0,9a)f(9)d0 where the 


^ *n* 

kernel k(9,9a) takes a different form depending on whether 9a 9r or 

a 

9a>S _ Qr. 

a 

A numerical solution ■*'0 this equation given a measured pQ(9r) has been 
implemented in the CSU CDC 6400 computer, PROGRAM FREDHOLM. Program listings 
are given in Volume II. 

The field proceudre for this technique involves taking a series of off- 
angle photos of a small plot (Figure 5 and 6). The probability of gap in 
each photo is determined by overlaying a transparent dot grid and recording 
the proportion of dots which do not intersect a foliage element. The vector 
of the probabilities of gap serves as input to a computer algorithm as discussed 
above which calculates the distribution of angles within the original plot 
(Figure 7). 

E. Diffraction Pattern Technique 

The diffraction pattern technique is much slower and more tedious than 
the Fredholm integral method, yet is still relatively easy and rapid when 
compared to either the point quadrat or orthogonal techniques. A limited 
localized distrubance of the canopy is encountered when field photographs 
are taken. Subsequent data reduction requires several photographic and 
measurement steps. 



reproducibility of the 
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Figure 5. The field procedure for the 
Fredholm technique involves making a series 
of off-angle photos of the canopy. 




O' Z 30 ' 60' 

F igure 6. A full set of multiple view angle photos 
consists of all 10 increments between 0 and 70 Z. 
The probability of gap in each photo is determined 
by use of a transparent grid. 




Figure 7. Distribution of angles for the 
entire canopy is calculated through the 
solution of a Fredholm Integral equation. 
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This technique requires a horizontal field photo of a portion of the 
plant canopy silhouetted against a white backdrop (Figure 8). High contrast 
film is used to enhance the contrast between the plants and the backdrop. 

The negative of the field photo acts as input to an optical diffractometer 
(Figures 9 and 10) which generates a unique diffraction pattern dependent 
on the angles of the foliage elements. A high contrast photograph is taken 
of the diffraction pattern, and the negative measured using a photo cell 
densitometer (Figure 11). A wedge blocking filter is used on the densitometer 
in order to determine distribution of angles in the diffraction pattern. 

This distribution, in turn, serves as input to a computer algorithm, PROGRAM 
PROP, which solves for the distribution of angles in the original field scene 
(Figure 12). Program listings are given in Appendix B. A mathematical 
convolution of several scene angle distributions then yields the overall 
distribution of angles for the area. 
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the Fourier method are taken 
of silhouetted portions of 
the canopy. 


Figure 9. The high contrast 
field photos act as input to 
a LASER diffractometer which 
generates a unique diffraction 
pattern for each scene. 




F igure 10. The input is placed 
between the LASER and a plano- 
convex lens. The quality of 
the diffraction pattern gener- 
ated is primarily dependent on 
the alignment and cleanliness 
of the optics. 



information in the diffraction 
pattern is achieved by rotating 
a narrow wedge blocking filter 
attached to a photometer through 
half of the pattern. Program 
PROr reduces this informatior to 
a density function. 
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Distr 


Diffraction Pattern 


Input Scene 


Figure 12. An example of Fourier analysis for determining the leaf 
angle distribution of a wheat scene. 
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III. DATA COMPILATION 
1. General Description 

The principle field data collected by TAMU/GSU for the canopy modeling 
effort consists of periodic canopy reflectance, intensive leaf area index 
(LAI) measures, extensive LAI estimates, individual leaf transmission measure- 
ments, and canc geometry photos. Relatively complete data sets are available 
for March 20, 1975 (Tillering Stage, TAMU), April 23, 1975 (Jointing Stage, 
TAMU/CSU), and May 20, 1975 (Heading Stage, TAMU). Less complete data sets 
were collected on November 24, 1974 (Winter tillering stage, TAMU/CSU) and 
June 26, 1975 (Ripening stage, TAMU/CSU). The follov'\)ing two tables summarize 
the data sets. A more detailed presentation of each of these data sets is 
included in the remaining four parts of this section. 


FINNEY COUNTY DATA SUMMARY 
(Collected by TAMU/CSU) 


I. 


March 20, 1975 


Time; 


LAI: 

Canopy geometry: 
Leaf transmission: 


Tillering Stage 


Field 416 

Plot 1 

Plot 2 

Plot 

1100 hrs. 

1045 

1030 

1145 

1130 

1115 

1300 

1245 

1230 

1400 

1345 

1330 

2.07 

Fredholm field photos 

4.06 

1.31 


10" LAI plots: Field 

367 

369 

370 

414 

421 

Plot 1 

1.82 

.70 

.50 

1.46 

.49 

Plot 2 

1.04 

.41 

1.60 

.39 

.35 
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II. 


III. 


IV. 


April 23, 1975 Jointing Stage 


Field 416 

- Canopy reflectance: 




Plot 1 

Plot 2 

Plot 3 

Time: 1000 hrs. 

To?T 


1115 

1130 

1145 


1200 

1315 

1345 


1400 

1715 

1730 


1800 

- LAI: 5-13 

5.36 


6.15 

- Canopy geometry: Fourier field photos 




- Leaf transmission: Not token 




- 10" LAI plots: Field 367 369 370 

414 

421 


Plot 1 2.22 2.15 1.48 

8.53 

4.54 


Plot 2 2.78 2.43 8.45 

9.17 

3.60 


May 20, 1975 Heading Stage 


Field 

416 

- Canopy reflectance: 




Plot 1 

Plot 

2 

Plot 3 

Time: 0945 hrs. 



1045 

1100 

1115 


1130 

1200 

1215 


1245 

1300 

1315 


1345 

- LAI: 4.11 

5.32 


6.04 

- Canopy geometry: Fourier field photos 




- Leaf transmission: Green, yellowing, dead 




- 10" LAI plots: Field 367 369 370 

414 

421 


Plot 1 3.82 1.83 3.12 

5.65 

1.80 


Plot 2 3.22 5.68 9.76 

7.64 

2.16 


June 26, 1975 Ripening Stage 


Field 

416 

- Canopy reflectance; 




Plot 1 

Plot 

1. 

Plot 3 

Time: 

1115 


1000 

-- 

1200 


1115 

— 

1245 


— 

— 

1330 


— 

- LAI: 1.79 

2.17 


2.04 

- Canopy geometry: Fourier field photos 




- Leaf transmission: Dead 




- 10" LAI plots: Field 367 369 370 

414 

421 


Plot 1 .78 1.15 .94 

1.79 

1.02 


Plot 2 .85 2.00 3.14 

2.63 

1.08 



In addition, the diffuse to direct irradiance ratio and soil reflectance were 
sampled periodically for each date. 
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LANDSAT RADIOMETER DATA 
FOR FINNEY COUNTY FIELD SITE, KANSAS 
(Averaged Over On-Off Row Set-Ups) 


March 20 

Band 4 

5 

6 

7 

PLOT 1 

1052 hrs. 

.076 

.078 

.235 

.315 

1136 

.065 

.072 

.217 

.316 

1304 

.062 

.059 

.209 

.299 

1355 

.067 

.070 

.210 

.299 

1635 
PLOT 2 

.067 

.075 

.240 

.329 

1044 

.075 

.071 

.260 

.340 

1127 

.074 

.070 

.244 

.333 

1252 

.066 

.070 

.237 

.326 

1342 

.067 

.069 

.231 

.322 

1515 
PLOT 3 

.074 

.086 

.242 

.335 

1033 

.071 

.083 

.298 

.475 

1114 

.055 

.075 

.300 

.443 

1236 

.055 

.050 

.250 

.383 

1331 

.060 

.055 

.258 

.390 

1623 

.068 

.073 

.288 

.430 


May 20 

Band 4 

5 

6 7 

PLOT 1 

0945 

.037 

.028 

.197 .320 

1100 

.038 

.029 

.201 .328 

1200 

.049 

.038 

.205 .293 

1307 

.048 

.049 

.206 .279 

PLOT 2 

1018 

.027 

.023 

.193 .290 

1119 

.029 

.027 

.161 ,264 

1220 

.041 

.033 

.210 .307 

1323 

.043 

.038 

.193 .302 

PLOT 3 
1042 

.039 

.026 

.236 .333 

1138 

.038 

.030 

.217 .376 

1245 

.040 

.036 

.232 .392 

1345 

,047 

.037 

.251 .347 


April 23 

Band 4 

5 

6 

7 

PLOT 1 

1009 hrs. 

.023 

.027 

.249 

.381 

1128 

.047 

.033 

.258 

.409 

1320 

.058 

.058 

.258 

.405 

1738 

.054 

.029 

.270 

.416 

PLOT 2 

1042 

.039 

.025 

.266 

.401 

1146 

.043 

.030 

.261 

.381 

1340 

.049 

.035 

.275 

.402 

1738 

.054 

.025 

.321 

.503 

PLOT 3 

1110 

.040 

.025 

.262 

.391 

1203 

.045 

.031 

.266 

.395 

1407 

.049 

.035 

.276 

.408 

1754 

.055 

.025 

.338 

.519 


June 26 

Band 4 

5 

6 7 

PLOT 1 

— 

— 

— 

— 

— 

— 

— 

— 


— 

— 

— 


— 

— 

— 

PLOT 2 

1110 

.083 

.098 

.153 — 

1154 

.087 

.116 

.180 - 

1240 

.092 

.124 

.166 — 

1325 

.092 

.117 

.158 — 

PLOT 3 

1006 

.070 

.134 

.155 - 

1119 

.078 

.113 

.152 - 
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2. Data Set Presentation 

Reflectance Data and Associated Parameters 


March 20, 1975 

Tillering Stage 


Field 416 

Crop type: 
Height: 

Satanta Wheat {10" drilT; 
8-9 cm. 

EW) 


Chlorotic: 

Weeds : 

Soil condition: 
Wind: 

Green foliage 
0 % 

Moist 

12-15 mph NW 

PLOT 1 

PLOT 2 

Leaf Area Index 


2.07 

4.60 

Dry Weight {2‘ X 
Number of Tillers 

2' Plot) 

(2' X 2' Plot) 

58.10 gm 

— 

Live 


803.00 

— 

Dead 


0.00 


Total 

Average Tillers/Plant 

803.00 

— 

Live 


8.64 

— 

Dead 


0.00 

— 

Total 

Average Leaf Area/Plant 

8.64 

— 

Green 


37.39 cm^ 


Yellow 


0.00 


Dead 


0.00 


Total 


37.39 



PLOT 3 
1.31 
39.90 

1372.00 

0.00 

1372.00 

11.00 

0.00 

11.00 

60.09 

0.00 

0.00 


60.09 


PROBABILITY 
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0 10 20 30 40 50 60 70 80 90 


LEAF ZENITH ANGLE (Deg.) 

ANGLE{DEG) 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 
85 90 

P(X) .031 .032 .030 .033 .034 .034 .034 .035 ,041 .050 .056 .065 
.073 .081 .085 .082 .078 .065 .061 


LEAF ANGLE DISTRIBUTION 
FOR MARCH 20, 1975 
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April 23 > 1975 
Crop type: 
Height: 
Chlorotic: 
Weeds; 

Soil condition: 
Wind: 


Jointing Stage 

Satanta Wheat (10" drill; EW) 
28-35 cm 
2-5% yellowing 
0 % 

Dry 

Calm 


Field 416 



PLOT 1 

PL01 2 

PLOT 3 

Leaf Area Index 

5.13 

5.36 

6.15 

Dry Weight (2‘ X 2' Plot) 

142.24 gm 

148.70 

207.63 

Number of Tillers (2' X V Plot) 

Live 

907.00 

1179.00 

931.00 

Dead 

52.00 

75.00 

33.00 

Total 

959.00 

1254.00 

964.00 

Average Tillers/Plant 

Live 

8.80 

10.80 

8.40 

Dead 

.60 

1.60 

.80 

Total 

9.40 

12.40 

9.20 

Average Leaf Area/Plant 

Green 

100.84 cm^ 

83.48 

89.70 

Yellow 

33.21 

42.50 

19.16 

Dead 

40.99 

32.02 

31.49 

Total 

175.04 

158.00 

140.35 



4 
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Field 416 


G. 


May 20, 197 5 
Crop type: 
Height: 
Chlorotic: 
Heeds : 


Heading Stage 

Satanta Wheat (10" drill; EW) 
72-89 cm 
13% yellowing 
0 % 


Soil condition: Dry 
Wind: 10-15 mph EW 



PLOT 1 

PLOT 2 

PLOT 3 

Leaf Area Index 

4.11 

5.32 

6.03 

Dry Weight (2' X 2' Plot) 

197.70 gm 

254.00 

287.00 

Number of Tillers (2' X 2' Plot) 

Live 

— 

— 

— 

Dead 

— 

— 

— 

Total 

— 

— 

— 

Average Tillers/Plant 

Live 

4.80 

5.60 

4.60 

Dead 

3.40 

1.40 

2.60 

Total 

8.20 

7.00 

7.20 

Average Leaf Area/Plant 

Green 

95.26 cin^ 

68.15 

151.27 

Yellow 

12.88 

5.09 

18.25 

Dead 

68.13 

28.97 

74.94 

Total 

176.27 

102.21 

244.46 



PROBABILITY 
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0 10 20 30 40 50 60 70 80 90 


LEAF ZENITH ANGLE (Deg.) 

ANGLE (DEG) 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

85 90 

P(X) .035 .027 .010 .007 .006 .006 .011 .014 .043 .061 .061 .083 

.098 .112 .120 .115 .101 .059 .033 


LEAF ANGLE DISTRIBUTION 
FOR MAY 20, 1975 



rspo aMn LOCaTlf. 
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Soil condition: Dry 

Wind: 20-25 mph SN 



PLOT 1 

PLOT 2 

PLOT 3 

Leaf Area Index 

1.79 

2.17 

2.04 

Dry Weight (2‘ X 2' Plot) 

181.40 gm 

219.70 

206.70 

Number of Tillers (2‘ X 2' Plot) 




Live 

443.00 

624.00 

515.00 

Dead 

71.00 

96.00 

89.00 

Total 

514.00 

720.00 

604.00 

Average Tillers/Plant 




Live 

5.20 

3.80 

5.00 

Dead 

2.80 

2.00 

3.80 

Total 

8.00 

5.80 

8.80 

Average Leaf Area/Plant 




Green 

0.00 

0.00 

0.00 

Yellow 

76.38 

54.80 

78.23 

Dead 

70.72 

42.59 

66.46 

Total 

147.10 

97.39 

144.69 

NOTE: The above measurements were conducted on 

June 18 and 

19, 1975. 

However, 


due to adverse weather conditions at this time, collection of radiometric 
data was delayed until June 26, 1975. Radiometric data was collected by 
NASA, LEG personnel. 



PROBABILITY 
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ANGLE(DEG) 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

85 90 

P(X) 0 .017 .009 .003 0 .001 .005 .011 .017 .026 .046 .075 .099 

.124 .147 .159 .143 .081 .040 


LEAF ANGLE DISTRIBUTION 
FOR JUNE 26, 1975 
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PROGRAM LISTINGS AND CONTROL CARDS 
SRVC 

Program Name: SRVC 

Subroutines Required: BLOCK DATA 

LAMBTN 

SUN 

ETHRES 

LANGLE 

NRM 

SE7Z 

UTIL 

COP 

PDENS 

PGAP 

OPTICAL 

Narrative: 

The CSU SRVC (Solar Radiation Vegetation Canopy) (Oliver and Smith, 1974) 
is a stochastic model for simulating the interaction of global radiation with a 
vegetation canopy to determine its apparent directional reflectance. Input 
requirements for the model are: 

1. Latitude of the target area 
Longitude of the target area 

2. Time, date and solar declination 

3. LAI (leaf area index) 

4. Leaf slope distribution 

5. Soil Reflectance 

6. Diffuse/Direct Radiation Ratio 

7. Leaf reflectance 
Leaf transmittance 

Outputs from the model are: 

1. Direction cosines of the sun 

2. PH IT — the probability of interaction of flux at specified 

source directions 

3. Predicted apparent directional reflectance values for discrete 
wavelengths. 
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A thorough doscn-ptlon of the model fs given in the report by Oliver and 
Smith (1974). An abbreviated discussion of the model and its application is 
also given in the Applied Optics article by Smith and Oliver (1974). 

A detailed flowchart of the SRVC main program is included and program 


listings 

are well -commented. 

Control 1 

Card Input: 


Card 1 

(8A10) 

Title 

Card 2 


(13) 

Day 


(14) 

Year 


(12) 

Hours 


(12) 

Minutes 


(F6.2) 

Latitude) 


( F7 . 2 ) 

Longitude) 


(F7.2) 

Declination 


(12) 

Bandwidth 

Card 3 


(12) 

(11) 

(15) 

Number of wavelengths 


Number of constituents 
Initialization variable 


(15) 

(15) 

Number of samples desired 


Number of trial 3 desired 

Card 4 


(no) 

Number of canopy layers 

Card 5 


(8F10.5) 

Threshold vector for downward flux 

Card 6 


(8F10.5) 

Threshold vector for upward flux 

Card 7 


' 


(no) 

Number of angles in leaf slope 

Card 8 


(no) 

Material type 

Card 9 


(8F10.5) 

Leaf angle distribution 

Card 10 


(2F10.5) 

S factor: leaf area index 

Cards 7-10 

- repeat for 

each canopy layer 
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Card 11 

(SF10.5) 

Wavelengths 

Card 12 

{no,7Aio) 

Number of input optimal vectors 
Description 

Card 13 

(8F10.5) 

Optical vector 

Cards 11-13 

- repeat for: 



measured canopy reflectance 
total study irradiance 
total diffuse irradiance 
soil reflectance 
leaf reflectance 
leaf transmittance 





Mathernatical Model 

The Monte Carlo model assumes that the canopy is composed 
of nonhomogeneous. layers of Lambertian surfaces of known optical 
properties, statistical composition," and geometric arrangement. 
The model is presented schematically in Figure 1. 


^ SUM 


SOUirE VECfOiSS 




N \ 


- ^ V ^ K S ^ 


>' / 


'tSO®r LAYERS 


; a '-' 


I 


LAMBERTIAN 

RESPONSE 


/I / . 
_ \ 


'1 


. . ^ ^ 9 I M 1 1 a I I u 

• / *' I OF FOLIA 

/ ^ / _ f 


/ STATISTICAL ENSEMBLE I ✓ N ' 

OF foliage elements 






background^ 


Figure 1. Schematic of a Plant Canopy Approximated by Stratified 
Foliage Layers Containing Statistical Ensembles of 
Lambertian Surfaces. 
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Global radiation which reaches the top of the vegetation canopy 
is composed of direct solar and diffuse sky radiation. The partitioning 
into fractions of total radiation is determined from the position of the 
sun, i.e., time of day, latitude, and solar declination, and from mea- 
sured global and diffuse sky radiation distributions incident on a hori- 
zontal sensor. The direct solar radiation component is treated as a 
point source whereas the diffuse fraction is divided into source sectors 
of the local hemisphere. These sectors are formed by partitioning the 
hemisphere into 10 degree inclination bands and further subdividing 
tbv^se bands to form 20 degree azimuthal sectors. The interaction with 
the canopy of each of these initial radiation sources is treated independ- 
ently. 

Diffuse flux resulting from the interaction of global radiation with 
a canopy element or with the background become new sources which may 
further interact with the canopy. The dov/nv/ard directed flux is combined 
with the appropriate hemispherical band of diffuse sky radiation. Upward 
directed flux is treated in a similar manner as diffuse sky radiation ex- 
cept the direction associated with each sector is the opposite from in- 
coming radiation from that sector. 

A frequency distribution of foliage inclination angle is determined 
for each layer from the geometric measurements of the canopy and is 
integrated using Simpson's rule to obtain a cumulative frequency di.stri- 
rjulion. This integral is normalized and partitioned into areas of equal 
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probc'ibility. The domains of each partition have equal probability of 
occurrence and may be sampled from a uniform distribution. 

The next step in the model is to calculate the interaction proba- 
bilities within each layer for both incoming and outgoing flux at each 
specified source direction. Several potential expressions are available 
in the literature cited (Idso and deWit, 1970; Pielou, 1969; a-nd Nilson, 
1971). The following expressions from Idso and deWit have been 
employed: 

p rt T L-A.I/ s 

P = Ll “ s g(0 ) sec 0 J 

where P is the probability of a gap 
o 

g(9 ) is the mean canopy projection in the direction 
^ of the source 

0 is the source zenith angle 
r 

LAI is the leaf area index for the canopy 

s is the leaf area index of statistically independent 
incremental canopy layers. 

it is usually adjusted by optimizing the reflection prediction for 
either one view angle or a set of wavelengths since it does not change 
with either wavelength or view angle. 

A given source finds a gap in the top layer of the canopy if a ran- 
dom number is smaller than P . The flux in this direction passes 

o 

through the top layer and reaches the next canopy layer. The absence 
of a gap necessitates the determination of material type with which 
contact has been made. This is accomplished by sampling from the 
di=»t ributlon of canopy con.stitucnt ► The orientation of the leaf is 
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determined by sampling from the inclination distribution and a 
uniform assimuthal distribution. These parameters determine the 
direction cosines of the leaf from which the angle between tne leaf 
and the source is determined. The optical properties of the Leai are 
then utilized to calculate the flux exiting the leaf in all directions. Each 
sector of the hemisphere on the reflecting side of the leaf receives flux 
for each v/avelength according to the equation: 


where 


2 2 

1=1^ psin(6j^g) {sin 0^ " sin 0^) 

Ts 

I is the source spectral flux 


o 


p is the material spectral reflectance 

0 is the angle between the leaf and the 
LS 



are the inclination angles defining 
spherical band 


source 
the hemi - 


The solid angle sectors receiving reflected and transmitted flux are 
defined in the same manner as for canopy flux sources only extended 
•Jo Include the entire sphere about the leaf. The leaf is not necessarily 
horizontal so a sector receiving reflected flux from the leaf is not 
necessarily directed upw-ard with respect to the local vertical. Hence, 
the direction cosines of the flux sector are rotated to the local vertical 
system and the flux pooled v/ith the flux in the appropriate source band. 
Transmitted flux is calculated and treated in the same manner as re- 



flected flux. 


v/l;ich passes through a gap or is reflected or transmitted 
downward from an upper layer of the canopy interacts with the next 
lower layer. Flux which reaches the soil surface is reflected into 
each of the upward directed source'bands. Upward directed flux from 
a lower layer of the canopy or from the background reaches the next 
higher layer and may interact with it. The upward directed flux from 
the top layer escapes the canopy. 

The interaction procedure continues until the level of flux in any 
source direction within any layer is below a threshold value. The flux 
exiting the canopy into each of the bands is separately accumulated. 

The ratio of the flux intercepted by a sensor placed within one of these 
bands to the global radiation intercepted by a vertical sensor with the 
same field of view gives the canopy apparent directional reflectance. 

Figure 2 shows the calculated mean response surface in the 
visible wavelength region for zenith view angles of 5 to 65 degrees. 

The non -Lambertian character of the canopy reflectance is evident and 
there is a general increase in canopy reflectance v/ith increasing 
zenith angle. This variation of reflectance with view angle is signifi- 
cant for the pattern recognition process. Sensor scan angle corrections 
may be required or this variation might be utilized as an additional 
characteristic feature. The distortion in the response surface with 
view angle indicates that methods employing channel ratio techniques 
for preprocessing (Krieglcr, 1971) or for specialized classification 
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api->roachas (Pearson, and Miller, 19?2) may require modification in 
some circumstances. 

The model was evaluated (Smith and Oliver, 197 2; Oliver and 
Smith, 1973) for a Blue grama canopy with a leaf area index of 6.5. 
Figure 3 shows a comparison of the model results, labeled SRVC 2- 
layer, with measured data for a vertical view angle. Agreement is 
good except for the chlorophyll absorption band. Recent evidence 
(Breece and Holmes, 1971) indicates that the foliage surfaces are non- 
Lambertian in this region which will necessitate a modification to the 
model assumptions for regions of strong absorption. Off-angle pre- 
dictions of the model are qualitatively correct but do not display the 
same precision as the vertical view case (Smith and Oliver, 1972). 
Prediction results using the original Kubelka-Munk model of Allen and 
Richardson (1968) are also shown for comparison. 


OEIGlNAi; PAGE IS 
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Fig. 


3 Comparison of calculated Kubelka-Munk differential 
eqSttion model and the stochastic SR7C model pre- 
dictions with measured canopy reflectance _or a 
vertical view angle. 
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ORIGINAI} PAGE IS 
OP POOR QUAIM 











sgT 

£)x, xa-^Pj xa.(pu 
Ti 9 sse^. 


0 ^ 0 , r,a)=‘ ^ 4 >u/:ce(xj kJ 

X=/jAfS<t-US 

K- h Hlam 


g/»tafJorJ / 5 £^X 
Szftfjsess AT i^Afsa / oi* 
DouMt otutjo xetftteHcs 


(CPM~ ETHKas' 


fJWY)NS 4 siif!,-l/ 


rsST KAD/ATJ!^ L£]fsLX 
/Af AU. OffiOfiuAXD Uiyegs 

AM/^/sr ra-zesaouD oaiuss 

fifiD X 6 ^D XubiCAragS 




ASPSa-CKi.) != X. S<pUHCsC^Spg,KL} rr X/a [t»X77}(XS<^b)J Sit} 


3 S<m^i 


tCL =* /j NLAf >7 


XDXn = “/ 


TCfTAU 

fiiun KSAxecrBD 
FKom A Uifn^STifif} 
zsFtecrsM. <?«o 
OBTSCTSO BY B 
Va/rr/CAi. SSAJso^ 


XL= O 


iXL=XL+/ 


I JS<fiR =0 


pOtUfJtoAR-p UiyeHS 


JS<J>ga 3 SCJCH 




/CAU. PSflP \ 

V fHT,NUVr, 20 IW 


DErEEtIJftJS Tfte PKOPOETJOW OP Ll&Ur (PItT) 
toHlCM IfJTEBfttfTS IM UfSVEK XU AMO THE 
TiPe OP MATERJ*!. fr»jrvpe) Bit SV th/j 
Semites. 


XPAXP=o 


XPHZPsi xrHxi>-t\ 


ORIGINAK PAGB IS 
OF POOR QUA3OT 


l - PHT 


/CAt-L LAHGIB 

Vra? 7 Ai.,PA 3 :i.) 


■’P 6 ts Tri'S pRoppBnoM oP <;ap 

DerBSHiioa kanpoH t^htaj. ^ Pwii.^ 
. SXMU, C 4 > 5 L, SIMP, C 4 >sp_, 

KtPj Y'.^p SLP> ‘“••••^ 

YMATie,j XHATIT 












TZS ^ Zs(jj‘s<pg,)'^ C^Sl, 

- ~?4S sxt^i. -h rss ^ xpmz= 4 


SZi. = Xs( J’JcSS'j XPHXP^ 
SYL = Ys( 7s4>/i_, xpfixp) 
SSrjL=r ssfj’s^^) 


Oiftec.T/afJ C4>sxf>iss 
ap soujsce ptMX 


/CflU,WMSTM\ 
^Cft,nrti»2j Js^^y . 
\F«TjIPSR) / 


contPttrss XUiPta /nAritiy: topic a catrrAiosia 
.THa pAtfitiior CP piax pspijEcreo f Tx^^tsmirrsp 
t^p Yig‘ TH£ aouKca 


jyjc^gs/ ,xfjCT s 1 1 


TJ^SijiK. - XTStfR-t I 


XX ACTS o 


xxpcr = xx/KT + f 


KL^ a 


HJLsKL-P t\ 


lamixxeiCTjiiu'i^ XLAe-7BCz:iAeri}iL)*icii(jjs4>}t-0 


j<L,wLfit£ y^^j 


^ZRCTiZ 


0—- 


XPPZsf O 


XPPX^XPPX-i-i 


2 </4 f nrs^n, xppx) 















K^FLSCTe'P FLUSC /J lioUluluJAXOj T/iii*JS- 
mirreb Ftu.v ts c4PtetAxo. captpura TrfS 
tAiaex zsi^s of tug So»xce £/xercr/oAj 
OF tfs /"4«y, 


R£FLS,aep puxx. ts uPt^A/tOj TRtttwnttTeD 
FlttX K Otx^tPA/iP. COF7PUT6 T^G tHOeX 
xsf»^ rpe sottace tuAecnon op 
me petfii. 


zs 4 >J = 


I ~~ ficc^s (zm.))i!cezTa l 4/.??? 


KL^O 


KLFKL.i -1 


KL ~0 


- HKCtfieY*X-Xt.i 

■p UifnB(z,iit,) 


UR(f}UC{i-Z-X.Uf Kt} = C/e(AHJir*Z-Xl, fu,) 

* UimsOtHtd 


0R(XLH, XSi>SfKt:i - DRCXlFt^ Zs4S, Xi.) 

+ LHmSOjKt.) 


ZS4-Sj Kt.) = DR Cxi. XS-i^Gj t({.) 

4 LflmCZfKl^ 


j g:NU3i^ 


~~~r^ 


IS 

r- 


PLc^y PASSi/^c- thsohOH 
Ft G-Fip tfJ CPHOP^ •••• 





> 

KL 

= 0 


KLsf fcL-H 


^,ZL-i'l,JF<S‘li,Xl.)F PR(sLPI,ye^PifXt.) 4 P& It DR(rU,ysel‘Z,HL'i 


/Siu-ssra \ 



< j M ' 

'eSw.sVAxes \ 


[ii-.‘Nuiy J-- 















U— 


/tsncHfAiC’ 

Safi, G/ZOUfJD 


^CfiLL ernAgs\^ 


T£sr tZA^fATlOM Levels M ALL 
• aat>Jf/MARP f-AYCAS ACAMST- 
r/fAestYOLD VAfMeg aa/o ser- 

XCipfi X/fP/CAT'OAS, 


J's<i> R9 o 


|£S0^ = JV^+ / 


[^^o(/aArY-/, jrs<^) : a 


!Cf>LL urx.t ^ 
XfHH0j3)j/}V£C 


Afa/£ HSAf/ soft, vscroe yaitp aybc. 


o»Vf I soft. VBCraa, 
Mpur ? 


PltJD APAfOoM 
So/i, i/ecra/f (AVee^ 


fif^sc/3) : t 


fCAlL MRM 
(cO,t,3),i{MU{/,(l^ 


filo^e VectoA Aif£e.” 
jAsra IZG 


fCALi. UTIL 
.RVSC , /?& 


JJ ~ ! 


XL=:Ul~ftV ■*• I 


HL- O 


KL^KLTl 
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un ( /; KL') - Oi JlTj *0.") + 0:1} XL^ fit) * Zs (js^x^ * XK (JJ"-!) 





C^iNstj^usT^)-^ 

T 


Set PLUX LS-\fBL /A/ oauffc^^so /cftu. SST^ 
Lnyex PiJiYi-i P<i'^ so»xe£ 

T'fsS/? r<S iSKC. •fST X&^)C 

/f/^ejfroA 719 aexo. 



vuiT-H, Jspu,-U 



33 




7SST X/3D)ffT>OtJ ISY&LS Jp All 
UPLt^AD LJfYSXS AaAlfJST ...... 

TffZesPOLO unties AAO SET 
XQ-<^ tt XtJLucATOaS. 



CALL £TO)!6s\ 

•j ffs<^ux,*\j 




xosr; 

= ■5" 1 


0 - 



Xl^ o 

1 f- 


J 

t 

1L=.3L-H 



XStfX= I 


,35 






'CALL PGflP 
^2L,a54>%MTY?fij^ 
\PrtTjHWYj IIO.R j 


OEIGaNAD PASS B 
OP POOE QUAUTB 


P&=;~ pfiT 









S6 
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KL-JiL-l'l 1 


Un»b(xxticr, AJi-P = •>i,iJi/yiQixj:ficTjKL) * XK(XXS^/i-i) 



xpux- C>\ 


XPftX-XPitX 


ZNL= ~iq£NL.(;r3's4'8,XPHr) 


flSPOSCTSO PLUt fs DOU}fittiSi^b^ 75?/JATJ- 

PiiTrsb PMX /s upiapdo. co/npurs ws 

WDEX XS^S Tf/S S^QX€iE 

Dxascrfo/J e>p 7-ps p£.ux. 


REPLSCXea piMX IS HPioP^O , rPA»SfHT7SP 
PLax ts fiapiuuAjib, Tpe /hbPK 

xsips ^P rue saupes oixicxio*^ oP 
7US PLU-X, 


xs4z- I <acc»5fe^)t■)■‘^‘y^)*^:garg (*%- Mj>s(zuC) » uk. | ■» 

[JO. » 10. i 


KL^O 


KU=i<L-H 


KL-Xt+l 


PS(«J.AY-i'2-3:i.j SS<fSjlCL) = Ox(MLfty+2-3Xj XS.(iSjAy5 

Ox(f7aft-tZ-Xt,,XS^SjfiLy-l>!i(NLtiyi-a.-Xt.fXs4>Sjlil^i-U\mB(l^KI^ ^ + UflMSCzjKi-'i 
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SKvC^NP^JTtOLTPUT.^'ILMPt«Pli^CH^ 0 SET 72 «T^PEl» 

»TflPFfc*PUTPUT«TAPec*IN(l-UT» 

r,.,. SCL/'W flftOI«TlON - VFeFTATTUM CANOPY REFtECTANCE MOOtL 
C,«,. THS pmoGHap CAICmLaIES THE fippAUEM PiPECTIONAl REFLfCTANGE OF a 
C.... VFfifc-rATTPK CAWOpy AS 6 FunCtICN OF CAnOPV 6£OHEtRYt tpAF Rt^LeC• 
C»... IrtNCE ANU Tf ANSHlSbICN* en|| PEfLFCTAnGE» ANU CANOPY TPRAOIANCE 
C.... For a OTVEn SoLaR POSITION, 

C.... P.E. CLTVFP AMD J.A. SMITH COLORADO STaTE UNIVERSITY JUNE. 197<« 


f 

r. CnMNCN RLOChS ANU HEFFNENC<iS 

c 

C LAPEL tKTERNAl. Rtf E PENCE. S 

r 

C C* HLf'Cy nATA, LAMBTN. SUfi ETHPESf LANGLEt NPM, SETZ. UTIL* 

C A^(; COD. 

r 

C G> LAHa^’N. PUENS* AhC» OPTICAL, 

r 

r C<.' LA^GLE. POPAS. AMI HGAH. 

r 

C G‘ KThRPS. «5tT2. ANP- lAMiaTN. 

r 

r O'- LAKGLt, 

C L" OPTICAi, . 

r 

C Gi '4T PGAH A^U LAH3TN 

r 

C // LA^Gl.£ ANO LAMBTf. 

r. 


CPyPOr. /r J /rAY.Y^AP»TTvE*GLAT»Gt0N6*CEGtRAK0i«»NLAM.TM£TSl tTHETS?. 
)NNATtt*TPA (4) .NCP.INIT.OUHI Il3> • 

?CF ClH.CFHTn.CEMTB.Ce^'lOP.rEiHI .c> ?P1 *C042 { l4> « 

■^SI^LA1 .rCS|.Ar,STNnec«CnSOFC.COCH.SlKZ»LOSZtSlNA7#COSA7tLXS<LYS.LZS 
erHvG\/P?/CANKM< !7» .SKYIM( 1 7) .0 IF I « U 7 ) . XM aT 1 R (1 7) 

) .uvAi t f n 7) ,xMATaB<ro ♦xMAi.^rn?) .xmat3bu n .xmatstut) .rrh 7) , 
r?*LAH( j7> .Srl'RCe < .THpT.a Um .ZEnIThUO) 

Cr.-'MOA/rA/f ANGLF (3i3) .FLA<3.3»10) .SLAT (3*3> •FLAT130) .PHIT 13*3* 1 «) 
CrHMO\/rf /rp <4* 1 17) .UR(4»10. U) • ThRESO <1 0» * lGflO<4t lo ) t IPOU (A.] {») 

1 . rHRE'itj n 0) 

c '■ kO N / r e / « I Nl » cOSL * s I ^P . GO 'P 

Cr:yn*ij4/t 1 /nATAin (7*9) .xmu(1*»*R) .0(17*17,9) *nvEC (R) 
erWH0iv/rM«T/RTP»3) ♦WLAYrOPMf 10) 

4»EM)LC 

CPHftlrx 4Vf:r ( 17) ,XK l9 ) . SXl *SvU.S7l * XLF ♦ Y^F *Zl> 

) - X*- (Ir .18) ,VS(lfi.iH) ,ZS (io) 

4 .ENPeH 

Intepnai. Arrays 

DIMFNbirN jOHIO(tt) *VEGT(17) ,SIG( 17) *V(17*17) *C0a<l7*l7) 

UlMENSlnK C0V(1,- • 17.17) *COVm(j 7.17) 

CTMEASinN vLFA (i9j *YLFA(lP) .0 m( 171 *QM1 07) *H£FER(l7) 

OIMFNSION ojT,lft.l7) ,BIT8AR(l0,1 /) ,RBAH(10|17) 

GlMEASlnN FUP).0P{9) 
peal L^SxLvS.lZS.TNCLF 
UTEGER DAv,yFAP«TH,TM*ZDE6 
MOrt. CALI. SETC(a.»CAnRMCi) ,enoLG> 

Gall SETG<r.,«*VEC(l) *EN0BH) 


Call •!iFTC(*>...JOi!»lDCl) .OP(0) ) 

r....BEsiHHtHaL cGMTrtPLS tttBCSaolNG PAGE BLANK NOT 


PA6B B 
OF POOR QUAUTK 



It- IS I = 0 
I‘=TCr = ( 

If'Ut = E 

(f ,1 (.3) IFli.F 

It- (c(;k (c) .^r .'•■.t s'^Ot- 

If ( Jl- I ST ) ri^Lt FtA’{-l«-l) 

• GcNtFflL «I»'I LA^ION CON«TPaTNtS 

FF An t iFTUe * 100) JOBIC.OAY* YfrAP.TH.TM»&LATteLONG,DEC*NpANOS* 
)^LA^/,^^/«T.T^IT.^SAMP*^^QI/^u 
IF (f-t)F <<=) ,,i.t STOP 

PFAO ( TF TLt . 1 ^?) NLA'' 

0'« = 3''/NP ANftS 

wi- TtE (f.poo ) JOuIn,OAY , YFAP. TH, TI'I.GLAT ,r-LCNP,neCtBANDWfNLAM.NMAT. 
1lNIT*f. FA^tP♦^rpIAL,^LAY 

RFaRI (FILE. 1 01) fHMFSC S-Pt An ( IF iLF « 1 0 1 ) ThRESU 
tvt-' T It 1 f- . irPi ) T3'-iF<;b . T t-RESL) 

.‘JACiif.-E f(:p u I ITAL I7 a nCN AMD CONVERSION 
NCOLPsKPANO':* 1 
M AyF’=^LAV•^l 
Fn v = <; , «rEnfP 
Fc vn?=Frv/?.o 
Sl^f.Ov = cx^ (FO'/O?) 

Call nAaSET(InIt) 

X I 1 = I H 

X r ? = T 

T^^'F = >n•►(''T?/6^ o 
Cil A I =.-i.a 1 AreuTN 
Gl OAfc^Gl UNO'-CFLiT-l 
Cc :=LtcocEmj{ 

ANOi- *CEl) IR 

.‘-Ln PC^'ITIOr pfWAMPicPc 
C'-LL son 

wFITtip.c??) ■> LYS • LZS 
Zc ( 1 ) = L 7 c 

.''CFCF IciFNTx: FOH UTFFLSF RAOtaTIO^ VECIUpS 
. sE' <=rp/MAM: area raTtc for all diffuse bands 

UhA3 = . 1 , 

S T X, A C = fl , 

Cl" ? 1 = 1 iNBANnS 
STNelsPiNA? 

ALSHA^sflLPt-A^ + bnNp.v 

sinapsstn (alpha?) 

Xx (I)=StNa?<»STNa?-Sina]asINa3 

2 DM ( r*1 > =STF FOV«STNFCv/XK ( 1 1 

hFnc.(o.2i''p) (Ak(T) .I = 1.NPAmOS) 
wt[Tt C".21F) (D*/l ( 1 ) , 1=2 .MSnUP) 

.fC'RcE DTFErTION InCi l.'.ATION AnC-lfS 
Tit al= h , 

f HF T A ( 1 ) = ( H aNiiW / 2, ) - e ANOw 
Cl" 3 t = l»NPAN0S 

Thfta(t*j)=thFTa(t) + bANOw 

3 TOTAL=IrTAL*(*K(I)/SlN (ThFT/>(i*1))) 

1 het A ( 1 ) =CFPl02_ACUS (COSZ) 

C('N5 = l ZS^^TpTAL 
C r I, r = 1 . 1 0 

ZFM 1 f- ( T ) =rFP105-TMPTfi ( I ) 

jvtrTt(N.??3) thfta 

.LilAtCriON CnsiNES OF azimuthal SECTORS IN THE DIFFUSE rAnDS 
CEG?n=2r .*cEOTR 


r>r>o jj no .jonono 


Cr fO JS0«e2«MSr,UP 

Zc ( jSt,P^ = ?T^J <TMrTA ( jSCR) ) 

P'-Isi .«CF[T>^ 
cr f-o iPHi = i.i8 

X? ( jSnw, IP hI) =C nS ( T mEI A ( jSOc) ) «COS (Phi ) 
v<=(jSLP.IPHT)aCnS(THtTfl ( J*^0o) ) (PHI ) 

Ph i = rr I *Ufe'^20 

C...•CA^rPV i;f0^fTHY. fAch canopy layer is copposed of one optical 

C. . . .^'ATERIuL WHICH May bF SPECIFIED AND UniOUE GEOPETPTCAL PROPERTIES. 
....C«xrPY t-eOHFTPIC HaHaNTTERS fOnSIST OF (jILEaF ANGLE FdEQUFNCY 

')ISTPi^iUTIO^ Ft'lMcTlON CENOTcn BY XLFA AND YLFA (?)|.EAF AREA INDEX 

. . . .t)Ei' 0 lEn RY flat aNR. (3)CAN0Py density denoted BY SLAI. XLFA (DEG) 
....ANi'! ylFa post hE SoErlFIED At aN ODD NUMBER (NANG) OF fVENLY ShaCEO 
....POINTS. FLAl IS N0N-NE(;ATT Yf AND SlAI haNGES BETWEEN o AND !• 

CFI F= m.^CFDTP 

went 

cr xb'- iL = l *Nl A Y 
RF AH iF TLP . 1 n?) NaNG 
RF AD ( IF tlE. 1 0?) I MAI 
P IP ( ii ) = IpAT 

RF Af ( IF ILP. ] 0) ) ( XLFA n ) . YLc-A ( I ) . IS] *NANG) 

fil-An(lFTLE.IOl) S| aI(Iu.IN’At),FlAI(IL.ImAT) 

. . . . tNTFC'Nf Tf aid ^0RY.ALI2F T^F | EaF ANGLE FHEGUENCY DISTRIBUTION 
. . . .FLf.'-l ION USING SlMPSfiNc RUI.E--THIS IS TEMPORARILY DENOtED by F. 
....N-i F(.iiA|.LY bPaCfO InIFRVALS of F aR£ THEN OFTERmInEO aND OENOTEC 
,...BY Fli (m POINTS). THF TaBLf FLA TS USEd FOr RANDOMLY SELECTING 
...•uEcF inclination AnGlEB, 
or 7(.5 t = 1,nAmG 

.'= XLF A ( 1 ) =XLFf ( T ) «CFUTR 
Ms ( (N-.NG-l )/?.)*] 

N/iNQLE (iL.TwAT)rM 

Call Tri R' (N . <LF f • YLFA «0M,F ) 

WwiTF (fc.?3'?> IF (I) » 1 = 1 »M) 

Cr 3J ■ TANG=1.M 
110 Fl A I 1l » THAT » I AN(. ) =UM ( I ANG) 

. . . .i.ccyal (ZF The (nrut leaf FREf'UFMCY cisthisltion funCtion to Obtain 
....A density FI'NC'Ion F which Ic fheCIFIEO at m points. 

FTDT= . 

CT' 3il T = ).^lA^lG 
in Firi = HTnT*vLFA (I) 

C<' .312 T=l.p 

•”2 F (I) = (Y| FA (?*I ) ,Y|.FA (2H»I*n ) /fToT 
Cr HE t=i,NAN6 
15 XlFA ( I ) SXLFA ( T ) *CFRTC 

write (A. 2 Jr) IL.ImATyNANG* (vlFA(1) .YLFA (I) .I=1.NAN6) 

WHIT£(A.?3n nAnGlE ( II . imat^ 

WRITE (6.23?) iFLA(Il.IMAT.Ii.I = 1*m) 

F = V - j 

wcttE (6.231) (F(I),1si,m) 

WP I IE (6.207) FLa I ( lL» IMAT) . SL A I ( I L » I M A I) 

. . . .rALCULA IF The mean hrclECTICm (OP) IN IHE DIRECTION OF ThE SOURCE 
...(THFta) of one UNTT leaf aRFA WITH INCLINATION INTLF. THE LEAVES 
....AT THiE angle ahf ASSLFED TO BE aZImuTHALlY ISOTrOPIC. 
cr 33‘ taNGLE-sI .NSOUR 
lNCl.E=-'=.*rBUTR 
cr 32 T*1,9 

Ir. DLF = InCLF*OF.Lf 

Call cop { Ih clf » thp T a ( I anglF) .np ( I ) ) 

r . . . . calculoTf yhf ^e^.m psCs;ecttlm johM) in The direction of the source 



r. . . . (ThET* ) OF nrE UnIT leaf ApEa AVE«aGE0 OvER the canopy leaf angle 
. .|)E^ c T T Y function F. 

Call copm (f .op. opm n angle n 

r . . . .Cfli,ri,tLftTF TmF PR08AdILITY OF A HIT (PHIf) FOR A LIGHT PAY WITH 
r....‘;Cloct OTHEcTION ImETA. 

Call PDF^5(lL.I^'AT»lA^Gu.E.OPM(IA^GLt) ) 

wp ITt (6.23's) CP. 0 pm f IaN'GLF ) .PhIT ( IL. Ima; , I ANGLE) 

3?(‘ COMINUF 
3='!t CONTlMJF 
WPIlE 

C. . . .pEFlEL r Ak.CE ANO THANSnISSION VFCTORS ape read for fach canopy 
r. . ..CCnfTI ILFNT, IN AnLJUlOM PEflECI ANCE VECTORS ARp READ FOR THE SOIL 
C. . . .pACKGRUUKI and TmF MFflFUREO gANOPY, THE MEAN VECTOR AND COVaRIANCE 
C....ANr cCPRFLATlOh waTHTCES APE CALCULATED aS WELL aS THE SQUaRE-ROOV 
C... .MATRIX wplct- IS SUhSELLENTi. Y USEO pCR MULTIVaRIaTE NORMAL 

c. . . c»-ui 1 TC wector Rampling. 

r 

c....*avflEngth'' ro bf sI'-^llateo 

PF AO I IFTLE . 1 01 ) ( XLAM ( I ) , T = i »N'LA'^) 

w Rile (6,2 0 \) (XL Ah ( I) . 1 = 1 .Nl«M) 

0 . . . .cc^ F r 1 1 ufnt opticai vEcroRR 

N('k = F'‘*AwaT*4 

ro 1 i L = i •NOP 

N. =(.-1 

IF(L.F.C.I) ^L=l 

lP(L.tG.2) call SplC(n..DATAlO(l*l ) •C(17.17.9) ) 

Call SFtc(p.,v(i .1 ) .V (17,17) ) 

CFAI UFtLE.10?) NVEC(Nl) .(UftTAll)(I,NL)»I = l*7) 
lF(L«ef:,1) NMsNvfcf'O) 

XwFO = NVFC (NL) 

wc[TE(E.2p?) (DaTaIij(I.NL) *t = i.3) .NVEC(N'L) 

I >'FC = NVFC <NL) 

C,...Pt'iD Ci-TTCAi CONST! Then r VFCtOpS 
Dr A 1 = 1 . 1 >/FC 

IF ( (L.NF.2) .Amu. (|_.NF .3) ) PFAQ ( IFILE * 101 ) (VECT ( J) » J=1 »NLAM) 

IR { (L.FO.2) .OP. (L.EO.3) ) REA0(IFlLE.ln3) ( l^ECT ( j) , J= 1 ,NLAM ) 
w-TTE (f. 251 ) (V‘^CT(.J) ,j=i .N|.AM) 

c. ..opitcal oonrtitufnt slns and Cross products 

CO a 

XM (^,M )=VMiJ( J,N|.) *vfcT ( J) 

Co 4 K = I , N L A 

V (^ ,K ) =V ( J.N) *VFCT (J) OVECT (K) 

IF (N ve:c (NL) .LF . 1 ) 60 TO 7 
C... .optical oONSTITUfnT CCVaRI«NcE MaTqIX 
Do c i 2 i ,tvLAM 
CO c ^=1 ,nlAM 

V ( I »w) = ( V ( T , J) - ( XnU ( I ,NL) *xmU ( J.NL) ) /X VEC) / ( XVEC-1 • ) 

IF ( I .'.-T . J) V ( I , J) = V ( V , I ) 

R CCNTInUF 

c ... .optical conrtitufnt n.ean vector 

Co k 1=1 .MAM 

6 XM. ( 1 ,M| ) SXNU ' 1 . NL ) /X VEC 

IF (I. .fcC.l) GO To 7 

Call matSCP ( V ,c ( 1 . 1 .NL) 'NLA'*) 

7 IF (L .NP . 1 ) GO To « 

C. .. .measured canopy RFflfctanCF wtCTQR FOR COMPARISON WITH MOOFL RESULT 
CALL oTtH vMJd .ml) .CANRM) 

Co 3fc • T = ) ,N lam 
CO 36' J=1.NLAM 



^ o o o ^ 


CCVi*-(I.j> * V(I.J) 

« rtc I Tt ( ^’•204) (XvU ( I .^L) ’ T = 1 
IP (WrC <M) .LE. 1 ) ‘sO TO II 
(^CIrt(f:.2l)5) (DiTfiIij(I.NL)»T = i.7) 

Cr 10 I = 1.^LMM 

10 W4ITE (<^.2P] ) (V( I. J) . jri ,t\iLaN') 

wr; ITc (0.21 1 ) 

C... .OPTICAL CONCTITUF'NT CCKREL"TtON wATRIX 
cr 12 I-1.^L4M 
Dr <; 1 *^T 4'^ 

S' crw ( I , J) = v ( I . J) / (COPT ( V ( I * II ) <»SOWT ( V (J* J) ) ) 
wPT T£ (6.2P1 ) (COR ( I .w ) . J = 1 *MLAM) 

12 CtntKUp 

11 CTMINUF 

r T T t ( * . 2 1 (■ ) 

F Ht(5 lODP 

ic j 0P = (I 

cr 7000 IPa.vP= 1 .iVSAviF 
Cr '■000 iToi m = 1 . M fp I A(_ 

C-* r Ti'^rZscECflMJ ( ApG) 

IP ( ( UPT Al. ,r,r . 1 ) . aNu- ( ISTOH.Eo.ri) ) fiO Tn 9050 
c . . . .''Er F RATE HAMjvJ^' nlPtCT AN(T OtFEUPE IMRAJIANCE VgrTORS pROM Th£ 
(■....^^PUT DIcTFIflUTIrMS A?’SU''''TN'3 THEY ARE vulTIVARIatE NORMAL. 

*>•- i Tt (6.29?) 

Of ?E I=l»? 

CclL uT tL ( yhU ( 1 . I ^ . avEC) 

IF (N \/Fr: n ) ,le. 1 1 ro tc ia 

C.'LL (C ( 1 » 1 » I ) .XmL ( 1 . 1 1 • A VFC) 

'4 rtP 1 1 1 (6.206) (UflTAln(^.I) . J=l,7) « (4VEC(J) *J=1•^'LA^') 

Gf TC ( 1 5* 1 6) • I 
r....TCT4L 6KY IPRAC,! a'MCE 
15 Call bTTL(/v£C.c'<YlM) 

C-f TC ?c 

r niPFuSt tRpaoiance 

’6 CoLL LTtL (AV tCtOlPlw) 

~S CfNTlMP 

r . . . .rCN PIJTE PHCPORTIOM OF IRRahUNCE 'VhICH is direct AMO PROPORTION 

r... .whirl- IS DiPFOsE. 

C(' 4 '1 J= 1 . r I. AM 

SrUPCE Mfj) = (SKVlM(j)-OIFl"(J))/(SKYlM(j)«LZS) 

Dr 4C I-s1,pRAmDS 

UC, SfljRCE (T*l . J) =OTPTM (^)*XK (n / (SKTJY { J) ASIK (THETa «I*1) ) > 
write (6.200) 
or 45 Irl.FFOuR 


43 

c . • « . 


• • • • 
• • • • 
^ • 


Q f)C 


wc r Tt (6.20?) ( SoORCP ( T . J) . J=1 .NLA»0 

►'CRIJL^IE first (TOD) CCWN (v.vfll l'^YER (OR) WITH TNCIOEnT DIRECT AND 
DIFFUSE I IGHT. nOwN C'wELL RaDTATION Flux (DR) IS TNOtxEO FROM 1 TO 
MLAY I.\ A UnwN GOImG SPGUFNCf. URwARO OwELL RADIATION h'LOX (UR) 

IS p.CEXFD FROw 1 TO NlAY+1 tN URwaRD GOING SEOOfNCE. TmAT IS FOR 

fcr uh, layfp 1 IS Ihe layer Immediately ascve (he background, tme 

FLL< lA LAYFP A'LaY*! is ThaT which ESCaReS THE CaNOPY aNQ TOGEThfR 
wlfh TrP IMrU'FM PUiX DFTFR|'"IN'EP ThE CAnOPY reflecTance. 

CCN r I auf 


Call 6ETC(n.«PEFER( \) .REEPRf l7) ) 
Call SETC ( a. ‘DR f 1 , 1 . j ) .UR (4. lO. 17) ) 
Call SFTC (n. I'~U 0 ( 1 . 1 ) , IGOU (4« 10) ) 


Of lOiM J=l.N50i|R 
D(' 10 O K=1 .NL A.' 


ir.fi3 DR ( 1 . J.K ) scPOPCf ( J.K) 





n o n t-) 


r.,..fi.ET FUju LEVEL I^OTCaTCRS (OoWN'WflRn) 

C >LL £ r^RE? (N lAv*nS0LR»-1 ) 

....fcrAL I^rInF^'T fl'Jx REFLEc^tn F^nM A Lambertian surfacf and de- 

....TECrFC ■Hy a vtPTrCflL SENSOR 
Di' '5--0 jSrR = l»NSnow 

or: i^i'O kl= 1 »nlam 

RFFF« (Ki ) =fFFFR (Kl) *SCUPCE (jSnM.M.) *5IN (THETA ( JSOfl) ) *9lNF0'/«SINF0V 
1A-1() COM INL'F 
1=''0 Cf’NTlMjF 

*|CTTE(r.2A?) (HfFfk(KL) ^KL^i «NLAM) 

Fast loop TDACE'i lI<?mt ATTEnLATIOn through CANOpY 

....fli'( Passing thhoogh layeps tn a do'ynwahd direction 

PnOil CO\Tli,l.lF 

Cr ?fc|.)0 II =] .-^iLa Y 
Cr ?5 "0 jsrP = i Soup 
C.... f;PPCF t-Li'X level InuICatOR 

IF nC'iji: ( II . jSfR> .Fij.O. ) GO TO 2500 
C...-Oln Ll'-HT «TCIK£ Li^AF 

Call hl'AP (il * JSnw.-i ,Ihit.mtypE) 

IF (ir I r.EL.c) Go TO 2200 
cr 21-0 IHHIP=1,1F 

r. . . .Olct£c I lOK CpSInES of SCURCF SECTOR (LVLH) 

S^L = *c ( jcpp, loHTPi 

SVL - rc ( JCOR. IFhtP' 

S7L = ZStjcOR' 

Call LAHbTh (lLtjSOH.MYPE.-1 ‘NSOUP) 

PI-'Ci CrNflAUc 

Gr TC PAon 

C....GAP ENCCUMF.REC.i tn jjOwKwAPO omTH 
R2''ii Dr 2E50 k|,=1»NLA'*' 

22 = 0 Cf ( TL*l . jSnp .kL) =(1R ( Iu + 1 , J:soR*KL) ♦DR(IL.JS0R»KL) 

2A*'iO Call SFtZ( tl* jSriR.-T) 
pci-i.', cpntIi'Up 

Call ETNHtS (NLAY.NPCLF.-n 
2eo0 CPNTIM'f 

C . . . .RACK^.ROO^D peached - PEFLECIs LAMBERTIAN 

Call E Threp (NLAy.nSolc.-t 1 
DO 3cAf) jCpMr 1 .k>OUH 
Ci'LL Ll TL (vmU( 1, 3) .aveC) 

IF (L vEC O' .LE. 1) GO 1C 3100 
Call (\R»'(C(1.1 »3) .XwLd ,3) .AVEC) 

Sion CALL oTtl ( AVtC»='’) 

Op 3‘* Ti JJ = ?.aSpuo 
IL = NL*Y ♦ 1 
CP 3*.P0 KLcI.NLam 

3 A.'.' Lp (1 . JJ.KL) =UR(i . JJ.Kl) *PG («-L) «0W ( IL * JSOR ♦ KL ) *Z5 ( JSOR ) *XK ( J J- 1 1 
C^LL StTZ tf LAv*i .jSOF.-l ) 

3fiPO CPNTIMJF 

Call fThREf (N|. Ay.nSoLh,*!! 

P....FLLX Pushing through layers tn an upward direction 
C o A(5 )0 IL = 1 .nLa y 
CP aSi'O jSnP = ?*ASnUP 

P....CHEPK FLiix level Indicator 

IF ( IGiJiJ ' !L. JShR) .FU.O) 60 Tn 4500 
P..,.(UC LIGHT STPIkE LfAf 

CALL pGaP ( IL. JSo-T.*1 , IhIT.MtYPE) 

IF ( ThI T .EO.O) GP TO 4200 
CP -il'iO IPHli^sl.lS 

COSINES OF SCiJRCE SECTOR (LVLM) 



S>L = xc ( jcofi. IpuTP) 

Svi. = Y<: ( j«:nK. IPHTP) 

S’l = 7«:(jc(iP) 

CflI.L L { 11 • jSOrt.MYPE**! ‘NSOUP) 

4l'i0 Cr\TlM;F 

G r TO 4 4 0 0 

C....G4P Fl'CCl'NTFPtP TN UP^^flRO PATH 

4^1)0 C*' A£i;0 KL^I’^La'^ 

Lt(TL + l , jSrq ,KU =HR ( lL + 1 » JSoP.Kl ) *UR ( IL , JS0R,KL ) 

44r'0 Call 5FtZ (II ) 

4=,0'' Cf■^TI^UF 

CALL '-.Ti-ipF<;(NLAY.NSOLP,*n 

Call tTHRE' (■ ni.ay.N)Solb.-i) 
call c fwME': (NLAv.rxtSOLc » + l ) 
r . . . .p£r ''CLF THPri.i(^H i.AYEyS LMTTl HlUX E'^HAUSteC 
Of c.jn.) II. = l«ALaY 
OF 5000 j5fa=?»^S0Up 
IP ( lOOIl ( II , JSOt:) ,i\F . .1 ) r?0 TO 2000 
^nn'; CFMjivIjF 

Cc 5ijf;l IL = 2 '''‘LaYP1 

OF ^O-Il J5 pp= 1 «r.5F|Op 

IF (TG'.in(IL,jSOR>, .ME.fn GO Tn ?0 0 0 

5 I C F \ T I MJ F 

C....FLl< tXHAUFTEH If all Fnu«CFc--COMHUTE RFFLtCTANCE FOR THIS TRIAL 
Cf 5i:00 JSFR = ?iaSdup 
OF 5F' 0 KL = 1 tA'LA'''. 

R T r ( w'sCFi-M.) ="H ,M AY* 1 , jcOR.KL) <*0'«1 (JSOR) / REFER (,xl) 

F?'-'’ RT Tmah ( j.fod,Kl)=RtThAP ( Jcr.K.KL) + RI I ( JSORiKL) 

Cr[A't = FirC0Ari(ARFl-CFTiyE 

WFT TE (fc.2a-,) iSa'-ip, IT = IAL.0t1mE 

CF 5J-0 JFrP=p fX'SnUP 

/f'EG= tOF-loojc.ot, 

= 3:-0 4PlTEt6.2H4) >DfG. (PIT ( JCO-f.^L) **VL=1 iNLAM) 

If ( r ivLfT ( aPG) .1, r .DTIf/E*2. ) tio TO 6100 

^'.?0 CFMIrwUf 

C. . . .TRIoLb CnPRi.fcTE fOP THIS SAMPLE POINT 
FTPl Als** Tpt AL 
C-F TC F-3 0 0 
f-l';,! FtriAL=tTPTAL 

4PTTE (6.265) tS.axp.I IRIAL 
''•2«-'i Cf ^3oO JSPPsP.ASnUP 
C(' 63"i) Kl=l,NLAW 

F3 :0 RT (FfAR ( JSOP.NL) = RT TbAP ( JF(TR.ixl) /F TRIAL 

k»c T TE (6 ISamp 

CF 44C0 JSFP=2,^50UP 
zr fg=1(!c-Io*JSOp 

f.4'-0 V«= I Tt (6 .2-J4 ) 70 fG. (PITBAR ( jFOP.Kp) «KL*1 .NLAM) 

Cr 66 r .0 jFnR = ?»xS0LN 
CF fib' O KL = 1 .N'LaM 

RPap ( jSfiR.kL) -Rh TP (jSCR.KL) *RITHAH ( JSOR.KL) 

Cf f-S' ij KLL = 1.Ni>m 

GFco crv (s,bOR.KL"NLL' =r.Oy (^SOR.Kl.KlL) *F1Thap ( jSqR.kiJ *RITrAR ( JSoR.KLL) 
Cf BCf-O ki.s1.mLAM 
R6«'f, F I TPAJ ( jSrp,KL) =i . 

IF ( ISTCp.FC. 1 ) GO To 71C0 
7000 Ci'NTIMjf 

F c APRs^ FAMp 
C-F TO 71 5(1 



72^13 


7 1 r 0 F« T 

....«LL SAVPI.E pnlMS t^bTI^ATEn 
71'=0 cc- 72'-0 jcnP=?»NSnUP 
G'' 72'!) KLa]*NLi'‘7 

WPdP (jSnR.KL)*-RPAW(jSCR»KU /FCA)-r’ 

C(' 7'»'''3 J*:OPs? tA SoUR 
7rFG=10'i-ln<>J?OR 
IF ( FS a^P .LF . 1 . > GO TC 7600 
DC 7n''n 

COVUSOP* I . J) = ICOV/ (JSCW* I • J) -FSA<'<P*P8A!^ ( JSOflf I ) <*rtBAR ( jS0« •- J ) ) 

'/{CCAl^R.]|.) 

7 4 0 0 S TG ( i ) =«;QP7 (COV { JGORt I . I ) ) 

CO 7f ! 0 1 = 1 *NL Af/ 

CO 75i u J = 1 .M.Ak/ 

Ci-H (i,j)=rov(jSc;R.l.v.)/ (SIG(I)^>SIG(J) ) 

VkPlTh (6.2=7) 7 Op G. (rear ( J?.Op*KL) •KL*! .MLAI^) 

I P ( F i t'^P .LF . 1 . ) GO TC 790(' 
aRTTL Ih.PPP) 

DC 7 7; (1 l=1.t^)LA^/ 

V«,cMt( 6 , 2 PP) (C 0 '/{J^C=,I.J) .Jsl.NLAN') 


7 5/10 


7 u ,• 0 
7f«lT 


77 >0 


7 h 0 0 
7'-;<' 0 


wpITE(R«2^i) 

CO 7ti;'0 1 = 1 ■ 

viPlTt (P.EPO) 

CONT IMJP 

IF ( IF ILP.EU.S) 0.0 ro POOO 
s r-~p 


NL Am 

(Cop ( I ,w) . j=\ »m-ap) 


riAja 

^..0 FruMAriqMr!/i4X.i3.7x,l4.7v.2l2»6X,F6.?*7X,F7,p,5*»K7.2*PX.T2./« 

i5x.y2.7y.Tl .7«.T5.'Jv»I5«flX*T5) 

T'l FrP^'Al(aFln.G) 

■ 2 for M 01 (T10.7AJO' 

J'-O f'?pmOT (oll’.'oTi.aSOLAR R AO I At I ON/ vFGE T A I T ON CANOPY REF|ECTANCF. 
H_o.//,f. 4 X.oTNPUT nA1A«.//.lx«8A10»/* 


MODE 


po 

4 « 

HI 

7« 

OHI 

WHI 

u« 


yeah •a.iA.o. Time «*2I2»* HOURS*./. 


„ULIAK OAV w»| 3 ,-. TC.. 1 . . 1 ..^ 

LATlTtlDE = *.F6.2.'0 OEGREPS. t,ONGlTUOE » *.F7,2.* OpGREESo./, 

POLAR DEOL I^'ATION = O.F6.-3X* DEGREES**/* 

rAAu hIOTH oF OtFpLSE vectors = *.F5.T«* DEGREES**/. 

\ljMujFP OF YcLFNGI H BAM3P SIMULATED o.I2./. 

OF CANO.PY CCNSTITUcriyTS *.ll*/* 

K CiC-n ^DO NO. T(j IMTIA, I7E HA^DOM SECUENCE * **!=*/* 

^P amp =* . I S. / , 

^TR1AI = *.15./. 

\iay = *.yi. 


p,-.y F('Pr/4y (h»owavELEnGTHS S I MJL ATtO* • / • * 0 * • E 7 . 4 * 1 6Ftf , 4 ) 

2 2 FOPMAi (h»O....**7'MO.*NUM0EP op VECTORS = **12) 
pij F^RmAF(hi *.F7.4.16FH.4) 

2 4 foomai (*o mf am* . / . 8 X . 1 of 1 ? . 4 ) 

p 5 Fr = MoT(H.o CoVaHIANCF mAtPIX«) 

p.,t, FrPMAT (hiORaaQom vpCTCR GEnEpatEH FROM IhE *.7A]ri./.(* *»10E12.4)) 
2 >7 FoPiyaT(*OLAl = ni.F4,2.4y.**S = *.F4.2) 

2'f^ Ff'RMfl T (*aDtFFuSF vector coefficients*./. 

IS (O K *) ./, (QFm.4) ) 

FORMA t (oOlPRADlaNcE SCURCE vlCTQRS*) 

2 1 F r R M A r ( 1 H I ) 

pil fprvamhm) CoRpEiATIOM m.TrIxa) 

p,2 FORMA I (*0I1F 1 = **pFp.4) 



7K0 FcAl^f^sTSAi^P 
.,..aLl “iAJ/PLE POINTS pSTI^ATEO 
71^0 C<” 72^0 J‘?nR=?»ASOUP 
Df' 72' !) KL = ] »MLiM 

72;j 3 PpAR ( jSOR^KL) pRpAh { JSCR»KD /E'^Ai'-rj 
DO JSOPs? SOllR 

ZrFG=lOS-ln*JS0ft 


IF (FSapp.LE* 1.) ‘’0 TC 760P 
DO 74'‘fi IsItNLA^ 

7 7'U) CCv{wS0p1*i! J) stcOV ( JSCW« I" Jl -’FSAMpttpBAR ( JSOR»I) »RBAR C jSORt J) ) 
1 / (PCiAt^o-1 . ) 

7^,0 0 STGU)=SQPT(COVtJS0RtI»Ii 1 
or 7=1 0 1=1 tNLA^ 
or 75! ii JsltNI.AK' 

7U.M; COR ( I . j) =crV (JS0R» I <w) / (SIG ( II <>Sl6 U) ) 

7e‘iO WBITt XOPG. (rear ( JS0P*KL) ’KL=1 »NL.AN') 

IP (FSavd.LF* 1 ■ 3 GO TC 7900 
WRITE ( A.PRP) 

DO 7?r() i=1,HlA^/ 

77;Hi WP[1E(F,2P03 (CO'3 ( JSCP » I • J) . J=1 •NLAM) 

w=ITE(h.29i 3 
DO 7E:'0 I = I»NLAw 

wCITE(ft.2flO) (COR { I ,w) « Jsl *’'ii-AW5 

79''0 Cr^TI^UF 

IF ( IF- ILP-EO.5) PQ TO POOO 


l-'O 


STOP 

f. data FOPi^ATS, 

FrByAr(BAIr./,4X.T3.7X,l4f7y2l2»6X«E6»2«TX»F7.2«5X»F7*2«8Xtl2«/i 

l5x.l2.7y,Tl»7'(,T5.9v,l5.fl)<«T5) 

I ^ I fcrwA I ( aFl 0 .5) 

^ - 2 F^a^/A I (Tin »7A1 O' 

?>-0 FORWAT («!■», a3<»«S0LAR RAD I At ION’/ '»FGE T A T TON CANOPY REF| ECTANCF MOD 
lL^*«//*e4X*'«'TNPUr oATA««//*1v*8A10»/* 

ULIAN’ day year time *»2I2»** h01JRS*»»/« 

latitude = «,F6.2.« degrees. lONGITUOE = *9F7,2,* DpGREEStt,/, 

POLAR DECL I''JAtION = «»FA.?*<» DEGREES'** «/ » 

AAND width oF utFfLSE VECTORS = ^iFS-T'* DEGREES**/. 

\ljf/HEP OF waVcLFNGIH BAA'OS SIMULATED <»»I2»/» 

MJMbEO OF CANOPY CCNSTITUfwTS *.Il»/» ^ . 

K D[C-tT odd no. To initialize random sequence = *.is»/» 

nSARP =•4,15,/, 


?•(> 
r '** 

O'** 

7« 

><■» 


u« MRlAI = *»I5./. 
so nlAy = *.Il« 


FORMAT (*>0 wavelengths S IMULATEd** . / . *0* * F 7 .4 » 1 6F8 . 4) 

2 '2 forma ( (*»0 ... .**341 0 .*NUMBEH OP VECTORS = *»I2) 

2^3 FrRMAf(o -tt.F /»4, lfiFB.43 
2 A FCOMAT (<*0 MFAN*»/,8X,10p12.43 

? 5 FoRMaHAo Covariance matrix**) 

^.,(3 FoRMaT (ttORANOoM VFCTCR GENEdATED FROM THE **7A]n«/»(* *»10E12.4)) 
^>7 FORMAT (**0LAl = *> . F4, 2 .AX.-aS = **.Fa.2) 

2-8 Format (aODTFFUSF VECTCR COEFEICIEN'TS*./. 

I9(tt N **) ./, (qFp.4) 3 
2 «c, FORMAT (*JOlPBAr?lANC£ SCURCF vEcTQRS*) 

21 ii Forma T (1 hll 

211 FPRMA!(«() Correlation m*trix*>i 

p,2 FORmAT(«0DM = **.RFP.4) 


2?1 (» 0 f HPESO = o.iOFF.i/* THRESU = *»10F8.4) 

FPRMP i (*pnTREcTTON CCFINE'^ nF SUN *«3F8.4) 

553 FrfiN^ri r (*0 ThETA =«.1oFP.4» 

?p 7 Fr-Rf.'A r ( ///<► 2t ( 1 H . ) . 2 X fOCftNOFv GEOMf T H Y* » 2 < ♦ pi; ( 1 H . ) // ) 

?.?H FrR^Aj (/* *,2S(1H.) ) 

210 FrCf/AT (« 0 L F AF AMiLt CCMPUTAtION? - IL c *,I1, 

1« l^aT = NANG = *»!<?./»* XLFAtYLFA** 

\/. (2X. I*F?.3) ) 

f-i\ Ff'RivA r («nNfl^GLE ( II. » T^AT ) r «..T2) 

PT2 F<''R«Ar(vO FLA =o.l0Ffl,3) 

?33 Fc«MAT(<»n F =*>,1.1F8.3) 

?55 FrRwAf(«0 nP =i».9F8.3*3x a ♦ , FP . 3 i 3 X , «PM I T = #,F8.3.) 

pE-i FCRMfl I ( bX . 1 OF 1 2. 4 > 

2P2 Fr^MA r (aopffer = i»*ftei3.4) 

P 03 Frui^A f (ftoRFFLf CTOnCF FOP SA '|HlE<* • T 3 » * f P I AL* » I 3 » SX » 

1 «trcp/Ru rxTlrA' TIX'E W(iS«,FS.l.'^ SECONDS.*) 
pJ 4 F'-Pi'^afC* 7 =4,lT.« rEGoOXo 0F7.3) 

r . . . .SAivPLF*. I3.4 OCNTaInS 0NLX«,I3»* trials.*) 

54h Fr-R^ a ( (*|i* . { lr( . W4 ypAN Rfi-lECTAncE FOP SAMPi.f:4 . 1 3 ) 

?u7 FCPVA7 (oOGRANP ^EaN FCP 7 =<>,13.4 DEGREES . * 1 3 ^ 1 1 OF 7 . 3 ) 
Fr:;N-aT(«p COVARTAh'CF MATRT*4) 

5 ►. 9 F r R A ) ( 7 X » 1 f F 1 2 . 1 

211 FrRi^ATlttO CoRfEi A1 ICN MAlwIXo) 

2 c 2 FrPvA r (1 X, ] 2” ( !>--) ) 

Et-P 


OOR QUALTTV 



-) O o 


s.i^oOUTTNP I-«"STN(1L.JS0H-"TVPE.I0IB.NS0UB' ,ppROPpIATE 

i E VE. EEOV ION .~0 OPtlE.U PBOPEBTIES 

■•■••;;::rt;r,/E 2 scae,epeo eu.«. 

Sx'l. SYL* -7L 
jSi-.q 

- LYS. L7^ 

- 1CI« 

nPtT •«-«•<) 
r uPti'^’K* 

p M T V Vi P 

r IL 

r NSrUP 

r XN aTl ♦ xf^AT 1 r 

r XKATPT 

C XIvAOH. XKAT3T 

r /Ea I T f 

Q O'JTPuT ViiHTARlP^S 
Q nP(I»>i**^' 

r UP ( I *o*a' 

Crii'^CA f'UPoUT) *'^L.*SyL»SZL»XLF»^L^» 

Cn»pOA/rl/nU«'l'’;'vL7;0“';f;f.,:”lTl!71.«"ATJ«(ni.«v»T2TlWI. 

Cry.>tvA/f /n ^ .MuaY.OH'^'I") 

Cry»vUN/r^AT/'UP(3) WTRP(?tl7) 

PPAL lX‘^’LYS»LZS 


cOVJdCE 


OTPrCTION COSINES OF NORMAL* 


1 r. 0 


5 

7 

9 

% • • • 


PPAL 

CaTA -In2/1.57079^3?7/ 

>^ET nlPE^TlnN COSINES CF 
Xi_ = 5X|. 

YI. = SYl 

IF (jSOP.NP,!) G"' ■'■0 100 
XI. =l Xs 
Yl=L'^S 

rtT/^k 

.OA.nO. LEAP nPlEKlAllCA. 

,a>E CPtTCAL "1 ‘ = 
rEi'.crp'rr'uMLENiiv-iL.i 

IFIILIh.H.. _ mtypE.THETAi *PHIL1 

call Mp'.ulCP L ight'stpikes. 

SEl Slop- OF lE/'P YPlCh Lie»-i 

Cri = XL*»XLE*YL*VL'=’*^L*'ZLF 

IF(0CT.lT.o,) IsinE=llUP 
Cf'SiSsAoS (po^ 7 ) 

IPnOIP.EG.l) 60 TO 5 

. "';.^:^;V.L'.■‘sOP.«L,.C 0 S, 7 E.^H,J 50 P,-IHETAL,/l.. 

^%?.^.^p'iL'‘sO-.KL>.COS(7E^,^H,.EOP.-tHE^AL./l.. 

“VcPMCAl PPOOEPIIEE FOP I.EaF IVPE 
• 5^70 (io.jn.SOl.PlVPE 
cc lE M =1 .nlaM 


and 


1SIDE»1 (TOP) ’ -) (BOTTOM) . 


1 0 


R (Kl.) =XMATir(KL) 

I 5 T (KL) =<MATi T (“•L) 

QC TC 

PO Cr 25 K| =1 ,^LAM 
R (Kt_) sXmAT?P (<L) 

■35 T (kl)=XmAT?T(KL) 

GP TO aA 

•^0 CP 35 Kl=lti'iLAM 
R (KL) =XmAT3R (KL) 

■35 T (KL ) sXwAT-jT (kL) 

«o CPMIMJP 

C....UPPATE DTFFU'SE SOURCE'S WITH <;CaTTEREC HAQIATiON FLUX 
DP >^0 jjSCf = 2.NS0IIH 

IF ( I S IDF.EO.-'' ) Call BFlU'' ( th^tAL . ZEM TH ( JJSOR) ,H. T»fl«NLAK^ tPTPP) 
IF ( TS lOF.EG. 1) call EFLOX (TmETAL'ZEMTH ( JjSOR) * W » P . T * K'LAM . PTPP ) 
LP S(' K| =1,HLAM 
IF I ILJR.EG.l ) GO TO 45 

DC ( TL *1 . JJcnW.Kl ) =L)o (IL'»I * JjSOR.KL) ♦PTHP (2*KL) 

Lc (KLA Y*2-TL. •.JJ5<''R«'<L) =UR ( Nl AY + 2- I L ♦ J JPOR »KL ) ♦PTRP ( 1 *kL) 

G'' 10 =ii 

45 C‘-(r..LAY*2-U« IJSORf KL ) =0P C'l AY*?- I L » J J^OR * KL ) *PTRP (2»KL) 
LC(IL*I.JJCPR.K|,)=U‘^(IL'*'1*«JJS0R.KL)*P1HR(1*KL) 

«=0 CPNTIMJC- 
fi- TL'Ra 
Em; 



o 


SL-=aOLTTNE PFi.Ux ( TA , T BP «H »H , T AM ,PTHP ) 

GTVFN The TBPaOIAmCE P of a leaf I^CLINED AT TA THIS PROGRAM 
C5TFR''IN’FS the flux BEFLECTrD AND TRANSMITTED InTO A fOURCE 

C Pt\C >vHnSE 7ENITH ANGI E IS tRP. 

CI^E^SI'^^ P'RPC>*17) .H(17)»d(17)*T(17) 

DATA HI /3. 14 159?6S4/,P I 0?/l .570796327/ 

Ft ( X * y ) sCOc (TA)«(SIM(X)**?- fIm(Y)«* 2) 
F;(x)=4cOS(-l/(TAN(rA)«TAN(v))) 

F3 (X « y» 7)=5.*SIN' (TA) *ciN<x) (nEL*.25* (SIN (2.«Y) -SIN(2.*Z) ) ) /PI 
CFl= . '87266463 

T i=tbp-pEL 
T?=TPH+nEL 

IF ( TA .LF.Pin2-T?) Go To 10 
IF (TA.GF.PI02-T1) GO TO 20 
GO TC 30 
C....CA':F 1 

1l' XFi=Fi (T2 iTi) 

CO 15 K| =1 ,ML 1M 
PTBP(l,KL)aB(KL)<»H(NL)«XFl 
PTRP(?tKL)=T(l<L)»H(KL)«XF1 
'5 CONTlNUF 
B P T U H i\ 

C....CACE 2 

->C XF 1=K 1 (T?,T1 ) 

IF (TA.lp. 1 .55331 r,U TC 21 
Pcp = H 10-5 
C-o TC 22 
^ 1 HPP = F> ( tBP) 

22 xp3 = F J (pPP.TI.T?) 

Cl' 2b Ki =1 .N'LAM 

PTBP ( 1 ,kL) =H (KL) * (R (KL) 4.T (Ki. ) ) *xF3* 

1 (P(KL)*H(KL) *PR?-T (kL)*H(!<Li<»PI*T (KL) *H (KL) «PRP) *XF1/bI 
H1QU(j,kL)=i-(KL)*(T(KL)'*’R<Kl))*XF3* 

1 (r(KL)*w(KL) »PRP-P (KL ) *H(KLi<»PI*R(KL)*H ( KL ) *PRP ) *XF 1 /P I 

35 CONTlivUP 
P r T Ij H ^ 

r . . . .CAPE 3 

3C TP = HU'2-TA 

XF 1 =F ) (TfitTl) 

CO 3b K|=1,NLAM 
PTBP(1.kL)=P(KU«h(KL)«Xf1 

36 FTRP (?,kL) =T (KL) (KL) »XF1 
IF (T 8 *T?.LF .3.106) GC TO 36 

ppp = h r03 
G( 10 37 

5h pcp = F; ( (T<3*T2)/3.) 

’7 /F 1=F I (t2,T8) 

CFL=((TdP*TA) /2. ) -.74176493 

XF 3 = F j (dBP ,t8. T?) 

Co 4() «| = I ,NLAM 

PTBP ( I ,tfL) sHTPP ( 1 .NL) *H (KL ) <♦ (P (KL) +T (kL) ) *»XF3* 

1 (p(kL)Ow(K|_) *prp-t 'KL)*H(KL)<‘pI*T (KL) *M (KL) *PRP) *XFl/pI 
PTBH (3.kL) =hTPP (2.6L) *H (KL) (T (KU) *R (KL) ) *XF3* 

1(T(kL)«h(KL) *oRo*-R (KL) «h (KL) «PI*R (KL ) «H (KL) *PRP) *XFl/pI 
4(' C^^TI^UP 
Ff TU^■l^ 

EF r 


I’- 



Si SPOUT TNE LAN'GlE t IL • YPE » tHET AL • FH IL ) 

OUTPUT '"pamtom^i Pa^^bcJ^c^ SINl* COSL. SINPf ANo COSP ARE ALSO 

seJecIeo/ '’-'^‘-Ectance and transmittance vectors ape also 


IN-UT 

IL 

NTVFF 

f ANfcuE 
OUTPUT 
TheTaL 
FPIL 

TLF » Y| F. ZLF 

'T^L^ rOSU‘ S1\P, cC®P 

'rVATiP, x^aTjT, xmaTjr, vMaTjT* waTDR, XMAT3T 


CCr/vON/ri/miMPdn *CERTD.nUM7 (3) tCEpPi 
C^«'PC^/^4/^ angle (3»?) ,FL A 0 . 3 . 10 ) . SLAT (3 

CfPMCix/re/ciNLtCOSL.SlNP.cOsP 
Cr'"/vCN nu^3 (?P) .XlF t YLF»7 lF 
DETEPNUF PA^OOM LF**^ CRIEmTaTion. 

Ff.=N/iNGuE (TlfMTvPp) 

XTsPAAF (0. ) 

A 1 = 1 . ♦ (pM-i . ) »X t 
I x = V I 


O) fFLAI (3t3) 


.PHIT(3»3.10) 


iFClX.fO. N ANGI E(IL»NTYPE)) TX=1 x-1 
IXP1=1X*1 

^CPNT inclination and aZimuTh* respectively, 

STNU = SIn (ThFTaL) 

C(’«L = CCc (TpeTaU 
STNP=S1n(PpILi 
Cr?P = rcc (PHIL) 

CCVPuTt I Eaf NOHMAL DIcECTIOm cosines 
XLF = -S lNL<*rnSP 
Y| =’ = -S INL*«INP 
2LF = 0C5| 


select RANOCt" LEaF PEFlECTaNcE 

Call optical (mtypd 

FRTUPn 

E^ n 


AND transmittance 


vectors. 




u c 


L 


SLHROuTTNE cop (ALPMA tSETAtOP) 
r 

C....THIC PflOfiRAM CALCUI.ATEC TH^ mEAN PPOJECriOK OF A UNIT l.EAF APEA IN 
r....Ti-F DiPErTION OK ThE SCURCF. THE LEAF IS INCLINED AT *N ANGLE 
C....ALP‘-A fiND IS ASShMEU TC 0E AtImuTHALLY ISOTROPIC. THE SOURCE 
C. . . .DicECT ION I« AT aN AZIKUTH Op ZERO AND AN INCLINATION oF brITA. 


r 


C 


CCMMON/rl/DUHl 03) * CEP 10? 

CP = CCS(aLPi-a)»SiN»BETA) 

IF (ALPHA. LF.BPTa) RETlRN 

T>-ETao I« ThF LEaF AZIpUTH AmOlE at which CP becomes negative and 
IS IN IHF FIRST rIJaDRANT. TmE FUNcTION OP IS SYMMETRIC AND HENCE 
IS AVEHArtED OVER LEAF aZIMUTw aNGlFS OF n TO PI PADIaNS* 

Theta r = ^^cOe(TAN(HPTA) /TAN' AlPHA) ) 

T/.NTUrTAN (TPETAr') 

Ce=cH*(i.*(TAMro-THETAo)/OEpI 02 ) 

Rp TI'fiN 
End 





■1 o n 


1 . 


SLfROUTTNF COOM (G,0P,CPM) 

r....This calculates the mean projection of a Unit leaf apea in 

^....THF DiRErTIrr' OF IhE SCURCF. rOPM) FOR THE SIMULATED CAnOPY. THE 
OP The CA^OPY ARE ASSiiMED TO BE AZIMUFHALLY ISOtROPIC. THE 

....OP function used In the calculation has been previously determined 
....FOR A Ul\/EN «=OURcE direction FOR LEAF INCLINATION ANGLpS OF 

c 5* i5* ...» 85 Degrees, g is the leaf inclination angle density 

^....^•L^CT1CN. 

r 

ClYENSlnN OP (R) ,G (V) 

CPy = o . 

Dr 1 1 = 1 tR 
1 CPv = CPH*OP { I ) <»G ( I ) 

RF turn 
£^ C 



onomo T ■> T’^nno 


it '.n-iMi 


SL‘9PQUTINE PDEN«:(IL,KTYPE*IaNGLE*OPM) 

.._Xi-is HPOPRAI^ COMPUTES THE oWoUfiBlLITY THAT LIGHT AT INCIDENT ANGLE 

theta(iakgle) interact^ with material type mtype withik canopy 

LAvpK !L. 

INPUT 
TL 

MyPE 

iangle 

rpH 
SLAl 
Pl.Al 

yheT » 
r OUTPUT 
C 
r 

CrwwCN/r2/rUM(3r-7) »THETAM0) 

COvwCn/oa/nANGLP CY*3) *FL* (3.3.10) «Sl*I (3*3) fFL-»I (3*3) «PHIT(3»3.10) 

acg=1.-(Sl Ainu. MTYPE) «OPM/<;lN (THpTA ( I AnGLE) ) ) 

IF (APG.LE.rt.) oO fO 1 
F'=AHe««(P'|Al(lL*MlYFF)/SLAT(U.. MTYPE)) 

Go TC ? 

1 P r = f , 

vtciTE (^. I'ln) IAnGlE 

Jon FOPvAr (IHO. o PO SFT TO 7tP0*»iI5) 

■3 CONTInUc 

phit ( iL . mtype. 1 aNRLF) = 1 .-PO 
P r T U ^ l\ 

EnH 


i» 



o n ■> n 


SL*5RClTTNE HGAPdL* I angle* IniPtIriITtMTYPE) 

C ThI« PPOrPAn determines if Am INTERACTION IS BEING MADE IN LAYER IL 

c uNCStrs THE material type of LAYEb IL. 

A 

c INPUT 
r IL 

c tangle 

r iriR 

r ^ L A Y 

MP 
H“T T 
OUTPUT 
IHIT 
-vTYPE 

CO"'MCN/r A/NAN gLF C?*3) .FlaO.S.IO) .SLAI ( 3 . 3 ) *FLAT (3*3) ,PHIT(3*3« 10) 
CCivvCN/rMAT/MTP ( 3 ) *NLAY 
IF (ILIR.LT.O) GO TO 10 
r. av^eh=nlay*i-ii 
GO TC 20 
IG IlayEh=TL 
?0 N T VPE=mtp ( ilAYEP) 

Ii-IT= . 

TfTCT = qAMF(n.) 

IF (PHlT (ILaYER’Y IyPe* IANGLE' .LT.TeST) GO TO 30 
lHlT=i 
TO PFTl'HN 

Er 0 


ORIGINAirPAQgi? 
OF POOR QUALOT 


^ o 


SLPPO^JTT^E ETwRfS (Nl^Y •NSO'Jo * 101*^) 

TI-IS PRCrRa^^ DfTERmInES (FOR EACH LAYER AND FOR ALL LIpHT SOURCE 
OIcECTlOvS) TF The SOURCE FLmX IS ABOVE THRESHOLD PEQUiRENENTS IN 
THF DlRErTIfiM INDICATED BV InIP. INDICATORS IGOD OR loOU ARE SET 
C ACrORL’INr-LY. 

r 

INPUT 

r M_ ■ Y 

r ^SOUu- 

r iciR 

C KL Ar' 


r 

CP 

r 

I-P 

r 

ThpE 

r 

OLTPUT 

r 

TGOD 

r 

/• 

IG('L 


CryNtCx/rl/DliH^ ( 7) .NLA^ 

cr.v^C\/r6/nR(At]0.l7) .UR(A*U'.17) .ThReSoHO) flGODlAtln) »IgOU( 4.]0) 
I .T^PESU (lAj 
C----(^Cl^^lv«HD FlUX 

IP (ILIR.GT.O) GO TO IP 

Nl ayER = N'Lay*1 

Dr ? In ,nlayfR 

D'' 2 wi=i »NccuP 

iGor; ( T « j) =n 

Dr 1 « = 1 ,N| AH 

IF (DRn.J*K) .LT.TmHESD(J) > OO TO 1 
I i'Gi I I * .1 ) = ] 

6r TO 2 
I CPivTlNUF 
? CCNTlM'F 
Rr T•.'R^ 

C UPhAPL F| UX 

’ ' COMIMif 

Dr 4 In ,Ni AY 
Df^ j = 9»NcpUR 

Ipnu ( i * J) =0 
Do 3 K-1 ,N'LAH 

IF (UR ( I . J*)0 .(.T.ThRFSL ( J) ) «0 TO 3 
IP01‘ < I * Jl =1 
Gr TO i* 

3 CPNTIMJF 

4 CCNTlAUF 

Rf IuRa 
END 


U C Cl O C C C C C 


SLapOUTTKE SETZ ( IL * I ANGLE . IniP ) 

r THIS PfiOGRA^ SETS THE FLUX (aND ITS APPROPRIATE INDICATORS) IN THE 

I' ICIR CIPfCTION At angle THFTrt ( I ANGLE) IN LAYER IL TO ZpRO. 

INPUT 
TL 

Tangle 

ICIR 
NLaN 
output 
DP 
L« 
iGOC 
IGCU 

CTHMCN/rl /PI.IM) (7) .Nl.AN 

CCNN0N/r6/nP (A * ) 0« 17) tURlAiiO.l?) *THRES(10) iIGOolAtlO) »IGOU(AtlO) 
IF(IClR.EQ.l) Gn TO In 

C DOnNinittRO FLUV 

Cr 1 k=i,M.AM 
1 DP (TL»IaNGI E»'^)=(>. 
igocmil.iangle) =0 
rftiirn 

c upward Fl LiX 

1 '' CPNTInUf 

DC 2 i<=l 'NLAM 
'd L‘R(IL«IaNGLE*O=0. 

IGCU ( IL, ^A^ GLF) -0 

fic’TURN 

E^n 





r — 

r 

r 

C 

r 

r 

r 

r 

r 

c 

r 


C-- 

1 '• 


1 1 
1 ? 


1 1 
c — 


21 

?5 


2 ^ 

C — 
3- 


31 

3? 


33 


SlPRGUTTNE optical (MTYPE) 

-THIS HHOpRAy SELECTS RANDOM l^-AF REFLECTANCE AND TrANSmITTaNCE 
VECTORS POP material type MTvPE. 

INPUT 

htyff 

^,VEC 

C 

OLTPUT 

xwaTiR, X^AT^T. XMaTjR. xMaT 2 T. XVATsHt XMAT 3 T 

CrMMCN/Ll/nATfil0(7»R),XMU(l7,R).C{17tl7,9)»NVEC(9) 

COMMON /C2/CANRM (17) .SKYlM(l7).DIFlM(17)fXMATlR(]7) 

1»*vAT 1T(17) tXMAT2R(l7) ,XwAT?T (17) .XmAT 3R(1T) »XMAT3T(17) «RGn7) . 
2XLAM(i7i , SOURCE (1 ft* ) 7) ,THETAnO) 

I=2«M type*? 

J=I*1 

GO-TG (1 0,20*30) .mTyPE 

■ -sEi.Ecr material type 1 vectors 

CONTINUF 

IF (NVtC ( I ) *LE. 1 ) RO TC 11 

Call N«M(cn«i»T) *xml (i*d »xmatir) 

GO TO 1? 

Call uttl(ymU(i.I) »xmatir) 

IF (N vtc ( j) .LE. 1) GO TC 13 

call NRm (C (1 *1 »J) *XML (1 ,J) * VMATIT) 

PFTURN 

Call uTtUxmU( 1 . J) «xmatit) 

RfTIJHN 

.-SELEff material type 2 v/ECTOoS 

COM INUF 

IF(AVEC(I) .LE.1) GO TO 21 

Call NRM(C(ltl*T) *X mL (1*I)ixmaT2R) 

G(' TO 2? 

Call UTTL (XMUd.I) *XMAT2R) 

IP (NVEC (J) .lE.1) go TC 23 

call NPM(C(1*1*J) .Xml (1* J) *xMAT2T) 

R T U R N 

Call UTtL <''mU(1 , J) «xMAT2T» 

PfTL'PN 

--select material TTpE 3 VECTOdS 
CCNT iNUF 

tf(NVEC(I) .LE.1) GO TC 31 

CALL NPM(C(l«l»T),XMt.(l*I) tvMAT3R) 

GO TC 35 

CALL UTTL(xMU(1,I) ,XMAT3P) 

lF(NVbO(J) .LE.1) go TC 33 

Call NRm (C ( i * 1 * J) ,XmiU (1*J)*yM4T3T) 

RETURa 

Call uttl (xmU( 1, J) *xmaT3T) 

RftURn 

End 


Sl.PoauTTNE MPV(C*7»X) 

C.....TI-IS program generates RANDaM SAMPLES FROM A GIVEN MUlTIVaRI ANGE 

C NORMAL distribution 

coMMON/ri/nUNMT) »nlam 
dimension y (17) ,Y (17) 
dimension C < 17* 1 ‘ ) 7) 

DO 10 jslth'LAM 
Y ( I ) sGAmSS ( 0. , 1 , ) 

1 CONTINUF 

Call VM(jLT (Y»C*X*NLAM ) 
call vAdD(V*Z.X,NlAm) 

HcTURN 

End 







SLRRCUTtNE MATSQR(V.CtN) 
DlwENSlnN V ( 17» 17) *C (17f 17) 

Dr 30 Ja 1 f ^ 

Co 30 Ial*^ 

C ( I .w) =0. 

Ip(I*LT,J) go Tr 95 

IF(J.NE.I) GO To 5 

C(If J)=v(I.J)/SoRT(V(l,l) ) 

GO TC 3o 

5 IF(I.^P;.J) GO To l5 
Sl,i^?=0 . 

IK=T-l 

K»1 

10 SL^'^ = SU^'2*C(I.t<)*C(IlK) 

K = K*1 

If(k.lf.ik) go to In 

IF (V 1 1 » j) -?UM2 .gE.O. ) GO TO 11 

CIF=V 1 1 . j) - cum2 

>NCITE(6.?m) T* J,V(I,^) ,SUM9,0IF 
?'-l FORMAT (*C*»2I1 0.3F1<». 6) 

C 

GO TC 3F. 

11 COMUUc 

C ( I tw) =RGPT ( V ( I . J) SL>2) 

GO TC 3n 
15 SL'^PP0=»<', 

I J= J- 1 
K= I 

“0 SLMPhOscL'MPRO^C (I.Kj^CUtK) 

K = K* 1 

IF(K.lE.IJ) 00 TO 2n 
C(I«w) = (V(T*J) -RUmPbC) /C ( J V)) 
cO TO 3o 
95 C(T.w)=P 
1C Continue 
RR UjPn 
End 




Block data 

CCWMON/rl/DUMOo' ♦CEC7R*CE«TD,CEMTB»CEPl02fCElPI»CE2PT ^ 
DATA CEnTE«CEHTD»CEMTP/.nl U53?93 , 57 , ?Si577'J5 » ,000P908p82l/ 
DATA CEPIC?,CE1PI,CE2PI/1 .5707963? , 3 , 1 4159265 , 6 .28318«;30/ 
EKD 



SLHROuTTNE T6LR(M, X. Yt XX. Z) 


c • • • 
r* • • • 
C . • • 

c • • • 

r 

c 

r 

c 

r 

c 

r 

r 

r 


.Ihic PHO^,PA^ FTNnS ThE INTFGoAl 7(X) OF THE J ‘ Jn 

,TC USIM6 SIMPSCNIS oUlE. THE INTEGRAL Z(X) 

. at yIzw-D twe table of z vepsus X IS then inverted to oeter- 
lilNE X a5 a function CF Z at M regularly spaced points along Z. 

INPUT^VARIABLES^^^^^^^ CF regularly spaced points along z 

X = SPECTFIED AT ?N-1 POImTR 
V = SPECIFIED AT ?K-1 POImTs 

TmE^TaBLE of X values Fop m regularly spaced points 
( p-1 INTERVALS) along Z. • 

7 s THE normalized integral OF Y AT X(l)t X(3), .... X(Zm-ij. 


OlvENSIrN »(lR).Y( 19 ),Z(ln).XIlio)*XX(lo) 

c. . . .SIMPSCNS rule integration 

1 ' Z ( I ) = 0 . •! 

cx = X(?) - X(1) 

?■! CG5t'J = ^*H 

wO = -«».J - ? 

3 " -1 = " 2 

A' l^S) = 7(J - I) * DX«(Y(J0) ♦A.«Y(J1) * Y(J2))/G.O 

tSr XTU) = XlJ?) 

X I n ) =x ( 1 ) 

r . . . .NCCM aL J 7.P I^TEGRAL Z(X) 


GA CO 70 J = 1 

7 r Z(j)= 7 fJ)/Z(M) 

C....FUC X AT N regularly spaced points along Z« 

XX( 1 )=X(]) 

FN = M - 1 


7 = 1 . 0 /EM 

JS = ? 

DC 12' y = 2*M 
ZT = K - I 
ZT = /T*F 

or. DC ] li' J =J?*M 

IF (/ ( J) - 7T) 11 0, 

loo G - t/T “ 7lJ " 11) 
XX (K) = XI (J - 1 ) ♦ 

c-r TO ns 

CONTINUF 

ITS = U 

i^A CGNTINUF 
RrTURn 
End 


100. 100 

/(7{.|) 
G* (XI ( J) 


- Z(vi 
- XI (J 


1 ) ) 
1 ) ) 


t* 



SlPRCuTtNE 

r THIS calculates The position of the sun 

c 

C INPUT 
C TTWE 

r GLA r 

r CEc 

'• output 
C cTkl<.T, coslat 

r cjnCec. cn«0FC 

c rrpf 

r ^0S7 

c cyNA/, cnsflz 

r LXSf lys. lzs 

r 

r TI’^E OF simulation (HOURS) 

C C'LflT IS «ITE geographical LATITUDE 

r GLONC Is SITf longitude 

C DEC TS coLAR DFCLiNATION 

c ►’IS sclar houp angle 

r COS7 IS COSINE OF SOLA® 7EN1TH ANGLE 

C C0SA7 IF COSINF OF SOLAR A71MUTH 

C L)^5* LYS» LZS aRE SOLAR niRECTiQN COSINES 

r 
C 

CO^>yC^/^I/^AY,YEAR*Tl^'E«GLRT*GLONG^DEC*DU^' (2A) t 
lCEDTR»CFRTn.CEMTR,UUy?nZ), 

?SINL A r . cCSLAI ,S l \DEC»COSnEC. COSH. SINZ ,COS2* SIN A7,COSA7,LXS.LTStLZS 
Real lXS*LvS»LZs 

|- = AhS(((12.-TlHF)«15.) *CE07o) 

SlNLAIrFlN (G lAT) 
crsLA I =ros (GLat > 

S INCitrxclN (CEc) 

Co5DEC=rOS (DEC) 

COSh=C0F (H) 

CoS7 = SlNLAT<»STN0tc + CCSLAT<»CnSDEC*CCSH 
STN7=3GrT ( 1 .-COc7«COSZ) 

CCSaZ= («InOEC-SINLAT«COSZ)/ (C0SLAT*SInZ) 

SINaZsSoRT ( 1 .-CoSaZoCCSAZ) 

LvF=SlN7*CnSAZ 
Ly5 = SIN7*SIN'A7 
L7S=C0S7 
RF turn 

End 







Sl 'r’ROl.TTNa VMIJLT ( XiYtZ*N) 
....vector multiplication 

CIMpNisION V ( 17) * Y ( 17» 1 7) (1 7 ) 

Dr 10 laltN 

Z (I)'". 

DC 10 J«ltK 

10 Z ( T) =2 ( T) ♦< ( J) «Y ( I* J) 

PPTUPN 

End 



SL'POOUT TKE VA00(X,Y.Z«N) 
DTwF^SIn^ V ( 17) .Y ( 17) «Z (17) 

cr 10 ixitN' 

10 Z ( T) = M T) *7 ( I) 

Be TLiSN 
Er.C 


FL^CTIOM GAUSSUfS) 

gfnemotf pakdom Samples foom the umvapiate normal distribution. 

X 1 sBAnF ( 0. ) 

XpaPANF ^n. ) 

C13CI^ (x,.2R3185*Xl)*SGPT(-2.*ALOG(X2) ) 
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2. COEFF 

Program Name: COEFF 

Narrative: 

This program calculates the alpha and beta coefficients of a linear cor- 
rection algorithm for sun angle. The form of Y = Ax + B results, where x is 
the measured response and Y is the spatially or temporally extended response. 


Control Card Input: 


Card 1 


Column 

1- 

5 

(15) 


Column 

6- 

10 

(15) 


Column 

11 

-80 

(7A10) 


Card 2, 

, 3 

, 4, 

etc. 


Column 

T 

10 

■TfTo.3) 

(LP) 

Column 

11 

-20 

(F10.3) 

(L) 

Column 

21 

-30 

(F10.3) 

(RFL) 

Col umn 

31 

-40 

(Flo. 3) 

(WAVED 

Column 

41 

-50 

(F10.3) 

(THZ) 


Number of wavelengths 
Number of solar zenith angles 
Title for computation 


Path radiance for sun angle THZ 
Target irradiance for sun angle THZ 
Canopy reflectance at sun angle THZ 
Wavelength band for reflectance 
measurement 
Solar zenith angle 


All possible non-repetitive pairwise combinations of sun angles for each 
wavelength are computed. 



T 


PRECEDING PAGE BLANK NOT FttilECI 
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3.0 DATAUSE 

Program Name: DATAUSE 

Subroutines Required: PLOTl 

IDIOT 


Narrative: 

This program is designed to reduce radiometric data obtained using the 
CSU J-16 spectroradiometer or other similar instruments. The program 
calculates reflectance from ratios of sample reading and reference panel 
readings. The program currently contains the option to plot reflectance 
versus wavelength for a given time period. The program is designed so that 
other options of data manipulation can easily be added in the form of 
another subroutine. The desired types of data manipulation can be selected 
and input at the start of the data section. 

Control Card Input: 

Card 1 

Column 1-3 (13) (Nl) Number of PLOT types 

Column 4-33 (1013) (N2) Numbers to select each data manipulation 

option in a specific order, i.e., (#l=plot 
of reflectance versus wavelength "PLOTl") 
position 1 cols. 4-6 will be first in 
sequence and so on. Up to 10 data man- 
ipulation subroutines can be created 
and inserted in the program. At present 
only PLOTl is included. 


Card 2 
Columns 

; 1-30 

(1013) 

(Ml-MlO) 

The number of times each data manipulation 

Card 3 

- until 

termination of data 

option will be repeated. 

Column 

1-6 

(16) 

(IDATE) 

Date in 6 integers 

Column 

7-11 

(15) 

(ITIME) 

Time in 4 integers right justified 

Col umn 

13-22 

(AlO) 

(ITAR) 

Target type 

Column 

23-25 

(13) 

(lANGL) 

View angle; 2 integers right justified 

Column 

26-30 

(15) 

(I WAVE) 

Wavelength; 4 integers right justified 

Column 

31-37 

(E7.0) 

(DIRRD) 

Direct radiation measure 

Column 

38-44 

(E7.0) 

(DIFRD) 

Diffuse radiation measure 

Column 

45-51 

(E7.0) 

(SMPRD) 

Target (sample) radiation measure 

Column 

66-79 

(A14) 

(ICOMA and 

Extra comments 



ICOMB) 




Termination of 
columns 1-6. 

Cards followin 
Columns 1-40 
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data in this format is indicated by inserting 999999 in 

I 999999 card (if option PLOTl is selected) 

4A10) (lTIT) Title of graph . . ^ ^ 

1 card for each graph desired under 

this option 
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4.0 PLANTS 

Program Name: PLANTS 

Subroutines Required: IMAGE 

PLOT 

STATS 

Functions Required: GAUSS 

Narrative: 

Program PLANTS is designed to generate an abstract plant canopy defined 
by a group of straight lines in three dimensional space. This configuration 
is represented by the projection of each plant element onto XZ and YZ planes. 
The elements in the scene can be controlled (IFIG1=1) in which the user 
supplies the X-Y coordinates of the plants position on the ground (XGRD.YGRD 
and ZGRD, where ZGRD=0.), the length of the plant (DIST), the inclination 
angle (THETA), and the azimuth angle (PHI). A second option (IFIG1=2) allows 
the user to designate plant position, as before, however, the inclination 
angie, plant length and azimuth angle are randomly dptermined. The 
inclination angle is determined from a normal disvribution with TMEAN and 
TSTD. The plant length is similarly calculated from DMEAN and DSTD. Output 
from either mode of operation includes a microfilm plot of the orthogonal 
views; a graymap printout of the same views; statistics of the elements 
in the scene; and a permanent file containing a 64 x 64 boolean representation 
of the two views. 

Control Card Input: 

Card 1 

Column 1-10 (no) (N) Number of plants in a scene 

Column 11-20 (110) (K) Number of scenes (cases) 

Column 21-30 (IlO) (KDATE) Date of the run 

Column 31-40 (110) (IFLGl) Flag for random/static modes 

Column 41-50 (110) (IFLG2) Flag for plant branching (unavailable 

at present) 
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Cards 2 and Following (IFLGl) - one card for each plant 
Column 1-10 (F10.3) (THETA) Inclination angle (deg) 

Column 11-20 (F10.3) (PHI) Azimuth angle (deg) 

Column 21-30 (F10.3) (DIST) Plant length (inches) 

Column 31-40 (F10.3) (XGRD) X coordinates of ground pt. 

Column 41-50 (F10.3) (YGRD) Y coordinates of ground pt. 

Column 51-60 (F10.3) (ZGRD) Z coordinates of ground pt.=0 


^ (IFLG2) - one card for each scene (K cards) 

Column 1-10 (F10.3) (TMEAN) Mean of inclination angle 

Column 11“20 (F10.3) (TSTD) Standard deviation of inclination angle 

Column 21-30 (F10.3) (DMEAN) Mean of Azimuth angle 

Column 31-40 (F10.3) (DSTD) Standard ..eviation of Azimuth angle 


3 + N cards (IFLG2) - one card for each plant 
Column 1-iC (F10.3) (XGRD) X coordinates of ground pt. 

Column 11-20 (F10.3) (YGRD) Y coordinates of ground pt. 

olumn 21-30 (F10.3) (ZGRD) Z coor.:linates of ground pt. Z=0 


Repeat cycle of Card 2 and Cards 3 + N cards format for each new sc?ne that 
IS input under IFLG2. 


Example Output : 

Figures 1 through 4 illustrate the current flexibility of Program Plants. 
Each of the following figures includes a printout of the line segments generated 
by the program and a print of the microfilm image of the computer-generated 
plants in both the XZ and YZ planes. The microfilm negative produced by the 
program then serves as an input to an optical diffractometer to produce the 
corresponding diffraction patterns. Figures 1 through 3 represent different 
cases of the static option (IFLG1=1) in which the coordinates of plant leaves 
are specified. Figure 4 is an example of the random case (IFLG1=2). In this 
example the zenith angles were sampled from a normal distribution and the 
azimuthal angles from a uniform distribution. 

In the examples presented here a starting x, y, z ground coordinate is 
specified for each plant leaf. The polar angles and line segment length are 
then either specified or randomly selected and permit the calculation of the 
end coordinates of the line segment. For example, in Figure 1 the first point 
is located at X=45, Y=50 (Z=0). Given a zenith angle of 45° , a azimuth 
angle of 270° and a length of 30 units yields an end set of coordinates of 
X=23.79, Y=50. , and Z=21.21. 
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Figure f. Output from Program Plants for Static Case 2 and optically generated 

Diffraction patterns. Note that within each view the border was masked 
thus, only those elements within the border affect the diffraction. 
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Figure 2. Output from Program Plants for Static Case 3 and diffraction patterns 
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Figure 3. Output from Program Plants for Static Case 5 and corresponding 
Diffraction Patterns. 
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Figure 4 


Output from Program Plants for Random Case 2 and corresponding 
optical diffraction patterns. 
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5. FRDHLM 

Program Name: FRDHLM 

Subroutines Required: KERNAL 

INVERT 

DENS 

COP 

MATMPY 

PCTAO 

SETC 


Narrative: 

A detailed discussion of the inversion of the vector of percent vegetation 


cover as a function of the view angle 

is given in Section II. 

Basically, FRDHLM inverts the following equation to solve for f(0), the 

leaf slope distribution: 
r 7T/2 

“ L 

K(0^,0j f(9) d0. 

r a 

Control Card Input: 



Card 1 

Column 1-10 (AlO) 

Column 11-20 (AlO) 

Column 21-30 (AlO) 

(LABT) 

(LABX) 

(LABY) 

Title of PLOT 
Label of X axis 
Label of Y axis 

Cards 2, 3, 4 

Column 1-80; Cards 2 & 3 

Column 1-30; Card 4 

(8F10.0) 

(G) Probability of Gap 

19 values using 10 columns for 
each value 

Cards 5, 6, 7 

Column 1-80; Cards 5 & 6 

(8F10.0)(RS0L)Input comparison 

Column 1-30; Card 7 


CURVE (can be blank) 

Cards 8 and following 
Column 1-10 (FIO.O) 

(GAMMA) 

Is a smoothing function ranging from 
1-10,000; any values for gamma 
may be selected within this range. 
One value per card. 
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6. PROP 

Program Name: PROP 

Subroutines Required: GRAPH 

Narrative: 

PROP accepts as input densitometer 
diffraction pattern of an orthogona 
distribution of leaf 

This program is referenced in Section II. 
Control Card Input: 


Card 1 


Column 1-4 

(14) 

(NRUN) 

Column 5-iO 

(16) 

(NDATE) 

Column 11-20 

(F10.5) 

(DERCT) 

Column 21-30 

(Flo. 5) 

(BADJ) 

Column 31-40 

(F10.5) 

DTEST) 

Column 41 

(ID 

(MTRAIL) 


Cards 2 & 3 

Column 1-80; Card 2 (16F5.1) (DATAE) 

Column 1-15; Card 3 


readings which wedge sample the 
view of a plant canopy. A plot of the 

original image is then generated. 


Number of runs 
Date in 6 INTEGERS 
Threshold value 
Base adjustment for aperture 
Minimum divergence test value 
Test for end of data (other than 
0 for end of data) 


Densitometer values 


slopes contained in the 


Repeat Card 1 and Card 2 and 3 formats for each successive group of 3 data 
cards until all desired data has been entered. End of data is indicated 
by a single card with some integer value other than 0 for Column 41. 
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